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Abstract. We study a moving boundary value problem consisting of a vis- 
cous incompressible fluid moving and interacting with a nonlinear elastic fluid 
shell. The fluid motion is governed by the Navier-Stokes equations, while the 
fluid shell is modeled by a bending energy which extremizes the Willmore func- 
tional and a membrane energy that extremizes the surface area of the shell. 
The fluid flow and shell deformation are coupled together by continuity of dis- 
placements and tractions (stresses) along the moving material interface. We 
prove existence and uniqueness of solutions in Sobolev spaces. 
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•^ ■ 1. Introduction 

r^ I 1.1. The problem statement and background. We are concerned with estab- 

C^ . lishing the existence and uniqueness of solutions to the time-dependent incom- 

pressible Navier-Stokes equations interacting with a nonlinear elastic fluid shell 
(bio- membrane) . Recently, there have been many experimental and analytic stud- 
ies on the configurations and deformations of elastic bio-membranes (see, for ex- 
ample, [5], [H], [Hj, gg, [H], [IS], and [27]), but the basic analysis of the coupled 

»y-s ' fluid-structure interaction remains open. The fundamental difficulties arise from 

the degenerate elliptic operators that arise as the shell tractions. As we detail 

CO ' below, the bending energy of the shell is the well-known Willmore function, the 

integral over the moving surface of the square of the mean curvature. The degen- 
erate elliptic operator arising from the first variation of this functional is a fourth 

^^ . order nonlinear operator that smoothes only in the direction which is normal to the 

moving domain. Our analysis will provide a well-posedness theorem and explain 
the interesting interaction between the shape of the shell and the flow of the fluid. 

2 ' Fluid-structure interaction problems involving moving material interfaces have 
f^ . been the focus of active research since the nineties. The first problem solved in this 

area was for the case of a rigid body moving in a viscous fluid (see [12], [20] and 
also the early works of [21] and [23] for a rigid body moving in a Stokes flow in the 
/\ ' full space) . The case of an elastic body moving in a viscous fluid was considerably 

3 . more challenging because of some apparent regularity incompatibilities between the 
parabolic fluid phase and the hyperbolic solid phase. The first existence results in 
this area were for regularized elasticity laws, such as in jl3] for a finite number 
of elastic modes, or in [2], [1], and [3] for hyperviscous elasticity laws, or in [22] 
in which a phase-field regularization "fattens" the sharp interface via a diffuse- 
interface model. 

The treatment of classical elasticity laws for the solid phase, without any regu- 
larizing term, was only considered recently in lOJ for the three-dimensional linear 
St. Venant-Kirchhoff constitutive law and in llj for quasilinear elastodynamics 
coupled to the Navier-Stokes equations. Some of the basic new ideas introduced in 
those works concerned a functional framework that scales in a hyperbolic fashion 
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(and is therefore driven by the sohd phase), the introduction of approximate prob- 
lems either penahzed with respect to the divergence-free constraint in the moving 
fluid domain, or smoothed by an appropriate parabohc artificial viscosity in the 
solid phase (chosen in an asymptotically convergent and consistent fashion), and 
the tracking of the motion of the interface by difference quotients techniques. 

In our companion paper ,6], we study the interaction of the Navier-Stokes equa- 
tions with an elastic solid shell. Herein, we treat the case of a fluid shell or bio- 
membrane. This is a moving boundary problem that models the motion of a viscous 
incompressible Newtonian fluid inside of a deformable elastic fluid structure. 

Let il C M'^ denote an open bounded domain with boundary F := dfl. For 
each t e (0,T], we wish to find the domain ri(t), a divergence-free velocity field 
u{t, •), a pressure function p(t, •) on r2(i), and a volume-preserving transformation 
77(t, •) : f^ ^ M^ such that 



nit)^Tjit,n), 




(1.1a) 


•qtit.x) = u{t,T]{t,x)), 




(1.1b) 


t + VuU - uAu = -Vp + / 


in n{t), 


(1.1c) 


div u = 


in n{t), 


(l.ld) 


v Def u — pld)n — isheii 


on r{t), 


(Lie) 


u{0,x) — uo{x) 


Wxen, 


(l.lf) 


r]{0,x) — X 


Wxen, 


(1-lg) 



where v is the kinematic viscosity, n(i, •) is the outward pointing unit normal to F(i), 
r(t) := dfl{t) denotes the boundary of fl{t), Def u is twice the rate of deformation 
tensor of u, given in coordinates by u*- -I- w'j, and isheii is the traction imparted 
onto the fluid by the elastic shell, which we describe next. 

We shall consider a thin elastic shell modeled by the nonlinear Saint Venant- 
Kirchhoff constitutive law. With e denoting the thickness of the shell, the hypere- 
lastic stored energy function has the asymptotic expansion 

E shell — ^Emem + £ Eben + 0{e ). 

The membrane energy satisfies 

Emem —J dS = "f timcs the surface area of F(i) (1.2) 

Jr{t) 

where 7 > is the surface tension, while the bending energy Eben is given by 



/ \{4^l + 2X)H^ -2fiK]dS, (1.3) 

Jrit) L -I 



'r(t) 

where H, K denote the mean and Gauss curvatures on F(i), respectively, and A/2 
and n/2 are the Lame constants (see, for example, [17). 
The traction vector 

* shell ^^mem ' ^ ^hen ~r ^\(^ } 

is computed from the first variation of the energy function Esheii] the traction 
vector associated to the membrane energy is well-known to be 

imem=lHn, (1.4) 
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while the traction associated to the bending energy has a simple intrinsic charac- 
terization given by 

iben = cT{AgH-2HK + 2H^)n, (1.5) 

where ct is a function of the Lame constants and Ag denotes the Laplacian with 
respect to the induced metric g on T(t): 

1 d 



^•^-^M-^^'"''^) 



1.2. Outline of the paper. In Section [2l in addition to the use of Lagrangian 
variables, we introduce a new coordinate system near the boundary (shell) and 
three new maps, rj^ , if , and /i, which facilitate the computation of the membrane 
and bending tractions tmcm and tbcn- A key observation is the symmetry relation 
(|2.7p which reduces the derivative count on the tangential reparameterization map 

The space of solutions is introduced in Section [3l and the main theorem is stated 
in Section ID Section [5] defines our notation, and Section [6] provides some useful 
lemmas. 

In Section[71 we introduce the linearized and regularized problems. The regular- 
ization requires smoothing certain variables as well as the introduction of a certain 
artificial viscosity term on the boundary of the fluid domain. Weak solutions of 
this linear problem are established via a penalization scheme which approximates 
the incompressibility of the fluid. 

In Section [51 we establish a regularity theory for our weak solution using energy 
estimates for the mollifled problem with constants that depend on the moUiflcation 
parameters. In Section [9l we improve these estimates so that the constants are 
independent of the artiflcial viscosity as well as other regularization parameters. 
This requires an elliptic estimate, arising from the boundary condition (jl.lb ). which 
provides additional regularity for the shape of the boundary. 

In Section (TUl the Tychonoff fixed-point theorem is used to prove the existence 
of solutions to the original nonlinear problem (|l.ip Uniqueness, following required 
compatibility conditions, is established in Sections |4] and Section [TOl 

The inclusion of the inertial term erju into the membrane traction tmcm will be 
studied in a future publication. 

2. Lagrangian formulation 

2.1. A new coordinate system near the shelL Consider the isometric immer- 
sion rjo ■ (r, go) -^ (M^, Id). Let B ^ T x (— e, e) where e is chosen sufficiently small 
so that the map 

B -.B-^R^ ■.{y,z)^y + zN{y) 

is itself an immersion, defining a tubular neighborhood of F in M"^. We can choose 
-ip-^ , a — 1,2 and -^ on ;B where ttt. 
lotes the normal derivative. 
Let G = B*{ld) denote the induced metric on B from K'^ so that 

G(y, z) = Gz{y) + dz (^ dz, 

where Gz is the metric on the surface F x {z}; note that Go = go- 



a coordinate system ^^, a = 1,2 and ^ on ;B where j^ denotes the tangential 
derivative and -S- denotes the normal derivative 
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Remark 1. By assumption, goa/s — a^ • a^; where ■ denotes the usual Cartesian 
inner-product on M" . Let C'ap denote the covariant components of the second fun- 
damental form of the base manifold T, so that Cap — — -/Vq • -Aj . Then, G^ is given 
by 

{Gz)af3 = igo)aP — '^■zCafi + ^ ffo Ca-^CfjS- 



Let ft, : r — > (— e, e) be a smooth height function and consider the graph of h in 
B, parameterized by : F ^ i3 : y i— s- (jj,h{y)). The tangent space to graph(/i), 
considered as a submanifold of B, is spanned at a point 0(x) by the vectors 



d 



dh d 



dy"' dy°' dy°' ' dy" dz' 



and the normal to graph(ft) is given by 



<y)^JH'iy){ 



G 



ap dh d 



d 



HV) Qya QyP Q^ 



(2.1) 



where J/i = (1 + h^aC^, M^pY/'^. The mean curvature H of graph(/i) is defined to 



Hv) 
be the trace of Vn where 



d 



(Vn)y = G(V^n, — -) for z,j = 1,2,3 



where 



= 75^ for a 
derivative. Using (|2.ip . 



dvji 



1, 2 and ^ 



dz 



, and V^ denotes the covariant 



{Wn)ap = Gfv\_ 






-id 



'' dz 



' dyP) 



{Gh)8p [JH^GTh,^),. + J^\-GYh,K, + r^3) 



(Vn) 



33 



G 



■K'Gl'h 



_d_ 



— J,j {-Gh ^,7^35 + ^33) , 



(^4 



J, 



'' dz 



d_ 
' dz 



-<7(5 



where Ff ■ denotes the Christoffel symbols with respect to the metric G. It follows 
that the curvature of graph(ft,) (in the divergence form) is 

H = -{J^'Gl'h,)j + JiM~Gth,T], + F^:3), 

or (in the quasilinear form) 



(2.2) 



H = -K'Gf 



(Jfl-y - Ju ^Gl hnh 



-iq/3 



h,ay + Gf^'^Fapiy, h, Vh), 



(2.3) 



where F^p denotes a smooth generic function of y, h and Vft. 



Remark 2. Note that Gh denotes the metric Gz=h{y)j o-i^d, not the metric on the 
submanifold graph{h). 

Remark 3. If the initial height function is zero, i.e., h{0) = 0, then H{0) = F^3(0), 
which is the mean curvature of the base manifold F as required. 
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Figure 1 . The maps 77^ and ij" 



2.2. Tangential reparameterization symmetry. Let Af denote the normal bun- 
dle to r, so that for each y G F, we have the Whitney sum K^ = TyT Afy. 

Given a signed height function ft, : F x [0, T) ^ M, for each t E [0, T), define the 
normal map 

77'':Fx[0,T)^F(i), {y,t) ^ y + hiy,t)N{y), Niy)e^y. 

Then, there exists a unique tangential map 77"^ : F x [0, T) ^ F (a diffeomorphism as 
long as h remains a graph) such that the diffeomorphism ri{t) has the decomposition 

77(., t) ^Tl%-.t)o rj^;t), r,{y, t) = r,^{y, t) + h{r^^{y, t), t)N{r^^{y, t)) . 

The tangent vector 77, „ to F(i) can be decomposed with respect to the Whitney 
sum as rj^aiy^t) = V%{y^t)'^T: + ^^,Kiil^iy,t),t)r]'^^-^ and hence the induced metric 
5q/3 = Via ■'n,i3 may be expressed as 



9ap 



{(Gh)Ka + h^i^h,a] orf 



Kc.r,p ■■-- 



•V 



Ko.r,p- 



(2.4) 



Note that Qi^a is the induced metric with respect to the normal map ry''. Further- 
more, we have the following useful relationship between the determinant of the two 
induced metrics: 



det(g) = det(Vo?7^)^ det(G,,) J| o 77^ = det(Vo77^)^ det(g) 
where Vq denotes the surface gradient. 
Remark 4. The identity \2.4^ can also he read as {iffYg — Q ■ 



o rj 



(2.5) 



Let y and y — (p{y) denote two different coordinate systems on F with associated 
metrics 

977* djf _ dr]^ drf 



9ap 



dy°' dyP ' 



9af3 



dy'^ dyP ' 



It follows that ip*g — g. Let H, H , K, K, n and n denote the mean curvature, 
Gauss curvature, and the unit normal vector computed with respect to y and y, 
respectively. Since H, K, and n depend only on the shape of F(i), these geometric 
quantities are invariant to tangential reparameterization; thus, the identity 



H ^ Hoip, 



K — K o ip J 



no ip. 



(2.6) 
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Similarly, computing the first variation of / H dS in our two coordinate sys- 

Jrit) 
tems yields 

(^AgH + HiH' - X))n] (y) = [(^AgH + H{H^ - i^))n] (y) V y = ^{y). 

By (I2.6p . we have the following important identity 

{y)= Ag(iJo^) {y) \/y = if{y) 



A^*gH 



and hence 



where by 



[Ag{HoTj--)]oTj- = AgH 



(2.7) 
(2.8) 



Ho,j-- = -J-^Gf 



5/37 



J^^G2'h,h 



pn-,s 



•lafi 



h^^ + G';^''F^p{y,h,Vh). (2.9) 



2.3. Bounds on 77'^. Let u'^ denote the tangential velocity defined by rj^ — u^ orf . 
Timc-diffcrcntiating the relation -q = rf'' ojf and using the definition of 77", we find 
that 

d 



oz 



^From the trace theorem, it follows that 

||"''||H2.5(r) < C"P(||/i||H3.5(r), WvWH^n)) I|w||h3(o) + 11/7*11^2. =(r) 



(2.10) 



(2.11) 



for some polynomial V. Since, 77^(2/, t) = y + L [u^ o ri'^){y, s)ds, it follows that 



||Vo77"(y,t)||ffi-5(r)<C^l + 
and hence by Gronwall's inequality. 



^l|//2-5(r)(l + ||Vo77^||Hi.5(r)j ds 



\\u'^\\m--(r)ds 



for t e [0, T] sufficiently small. Furthermore, we also have 



(2.12) 



(2.13) 



2.4. An expression for tben in terms of h and rj'^. Now we can compute ibe 
in terms of h and rf: the highest order term of AgH is 



{vi^a?T[v/^^^^'^^^(^^"^(<-^'r^^rGr/..^.)/.,.,)]} 



^/M{g)dy-' 

Since Gaf3 — {Gh)ap + h^ah^p, the inverse of Q^s is 



o 77 



det(e^) 



22 



(Gh) 



-(G/i)i2 - ft.,l/l,2 

(G,)n + /72i 



which can also be written as 



T"/3 



A" 



ga/3 _ (_i)-+'Xdet(G,0-i(l - <5„,,)(1 - V)/^.k7i,. 



Therefore, the highest order term of AgH can be written as 

1 



Vdet(5o) 



ydit(^A"'3^^/7.„^ 



,7(5 



O 7/ 
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where 



j^aP^S _ J- 



Gl 



{~\Y+^ det(G^)-i(l - 6^^){\ - 5^a)h.^.h^^ 



(2.14) 



-2 /^f^K^Sa 



is a fourth-rank tensor. 



2.5. Lagrangian formulation of the problem. Let r]{t, x) — x + L u{s, x)ds 
denote the Lagrangian particle placement field, a volume-preserving embedding of 
a. onto ri(t) C K^, and denote the cofactor matrix of Wri{x,t) by 

a{x,t)^[VriixA)]-K (2.15) 

Let V = u o rj denote the Lagrangian or material velocity field, q — p o -q the 
Lagrangian pressure function, and F = f o -q the forcing function in the material 
frame. In the following discussion, we also set e = 1. Then the system (|l.ip can be 
reformulated as 



m = V 
-,3 n ('„.^*^ — I'^fe, 



vl-iyialD„iv)l)^,^-{afq)^k + F^ 



a,v , 







{vD,^{v)\- qSDa^Nj = (jQ 



L{h)B4-Gfh,^,l) 

ht = B,{{-Gfhc.,l)) ■ {v o r,-^) 
V = uo 

77 = Id 



V 



in (0,T) X n, 


(2.16a) 


in (0,T) X n, 


(2.16b) 


in (0,T) X n, 


(2.16c) 


on (0,T) X F, 


(2.16d) 


on (0,T) X r. 


(2.16e) 



on {t = 0}x n, (2.16f) 
on {t^O}x r, (2.16g) 
on {t^O}x n, (2.16h) 



where Dn{v)\ := {a^v\ + afv f.), N denotes the outward-pointing unit normal to 
r, O is defined in Remark [5l and _B* is the push- forward of B defined as 

B.(7'(0)) = (i3o7)'(0) V7(t)cr. 

L{h) is the representation of isheii ■ n using the height function h. It is defined as 
follows 

1 



L{h) = 



v/dct(.go 
hL2{y,h,Dh,D^h) 



y^^ti^A'^^'^'h^^f, 



,7(5 



L'^''^{y,h,Dh,D'h)h^p^ 



where Li and L2 are polynomials of their variables with Li(t/, 0) — 0, go is the 
metric tensor on F. Note that imem is included in L2 since it is a second order 
operator of h. 

Remark 5. For a point ri{y,t) e r(i), there are two ways of defining the unit 
normal n to T(t): 

1. Let n = y^^ a"'" N where N is the unit normal to F. 

-i(-Gf /i,„^-f ^)] oT^- (denoted by [Ji^\~V:Al)] ° 



J,: 



2. Let n 

if). 
The junction is defined by 



e(-Vo/io?7^,l) =a^N. 
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Equating the modulus of both sides, by i2. 5]) we must have 



e = Vdet(.9)[(J,7i) o 77-] = det(Vo77")v/dct(Gh)o7y-. 
Remark 6. An equivalent form of \2.16e\) is given by 

ht^~ h^a{v O ri^'^)a + {VO T]^'^)z. 

This equation states that the shape of the boundary moves with the normal velocity 
of the fluid. 

Remark 7. For many of the nonlinear estimates that appear later, it is important 
that L{h) is linear in the third derivative h^affj- 



Remark 8. Without using the symmetry 112. 8\) . we can still compute AgH in terms 
of h and rj"^ by using |^.^[ j and i2.5\) : however, Li would then depend on Vg?/'' and 
thus lose one derivative of regularity, preventing the closure of our energy estimate. 

3. Notation and conventions 
For T > 0, we set 
V\T) = {ve L'{Q,T-H\n)) I vt e L\Q,T-H\n)')y 

V2(r) = |weL2(0,r;ij2(f7)) w^ e l2(0,T;L2(O))}; 

V^{T) = {v^L^{Q,T-H''{^)) vtEL^{0,T;H^-^{n))\ for fc > 3 ; 

n{T) = \^heL^{0, T; ij5-^(r)) ht e L^iO, T; il^-^(r)), hu £ L'^{0, T; H°-^{T))^ 

with norms 

2 



v 



Vi(T) 



II^IIl2(0,T://i(O)) 



I II l|2 
+ ll^t|lL2(0,T;ffi(O)')' 

ll"llv2(T) = II^IIl2(o,T;//2(o)) + ll«t|lL2(0,T;L2(n)); 
2 I IL. I|2 



ll^ll-H(T) 



L2(0,T;ff'»(O)) 



+ ll"*llL2(0,T;ff'»-2(Sl)) for fc > 3 



"llL2(0,T;ff=-5(r)) 



''■*llL2(0,T;_H2.5(r)) + ll''-tt|lL2(o,T;ifO-5(r))- 



We then introduce the space (of "divergence free" vector fields) 
Vy^^weH^ (n) a\ {t)wij = V i e [0, T] } 
and 

V^iT)^{w£L\0,T;H\n)) | a^(t)7«^ = V i G [0,T]}, 

where the cofactor matrix a is defined by (|2.15p . We use Xt to denote the space 
V^{T) X n{T) with norm 

ll(«»llx^-|l«llv3(T) + ll^llw(T) 

and use Yt, a subspace of Xt, to denote the space 

YT^{{v,h) eV^{T) xn{T) ht£L'^{0,T;H\r))^ 
with norm 

ll("'^)llyr "^ II("'^)IIxt + II^IIl°°(0,T;H2(O)) + \\HL=°(0,T;H'^{r)) 
+ ll'*t|lL°^(0,T;_ff2(r))- 

We wiU solve (I2.16P by a fixed-point method in an appropriate subset of Yt. 
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4. The main theorem 

Before stating the main theorem, we define the foUowing quantities. Let qq be 
defined by 

A<7o = -Vuo: (VMo)^ + i^[4^^K)^],fc»(0) + divF(0) in n, (4.1a) 

qo = i^iBci uo-N)-N- cri(O) on T (4.1b) 

and 

Ui^vAuo-Vqa + FiO). (4.2) 

We also define the projection operator Pij{x) : R'^ — ^ T!,^(j, f)r(t) by 



V^,{x) = [% - (J-2 o T^^)a'ya'jNk{x)N,{x)] = 6,, 



4Nk{x) a]N,{x) 



\atNu{x)\ \a\N,{x)\ 



Theorem 4.1. Let v > Q, a > he given, and 

F € L2(0, T; H'^{n)),Ft G L^{0, T; L^{n)), F{0) G H\n). 

Suppose that the shell traction satisfies the compatibility condition 

[Beiuo-N]tan^O- (4.3) 

There exists T > depending on uq and F such that there exists a solution {v,h) G 
Yt of problem i2.16\) . Moreover, if uq G H^-^{fl) C] iJ^-^(r) and the associated ui, 
qq also satisfy the compatibility condition 



CP: 



i/(Def-uo)- -qoSi 



Ni 



+ v{5^i - NiNe) (Deiui)l ~ ((VuoVuo) + (VuoVuq)'^ 
- {Su - N,Ni) [iy(Bei uo)i - q^b^ <,JVfe = 
then the solution (u, /i) G It is unique. 

5. A BOUNDED CONVEX CLOSED SET OF Yt 



N, 



(4.4) 



Definition 5.1. Given M > 0. Let Ct{M) denote the subset ofYx consisting of 
elements of {v, h) in Yt such that 

\\{v,h)\?Y^<M (5.1) 

and such that v{Q) = uq, h{Q) = and ht{0) = (-Bo)*((0, 1)) • Uq- 

Remark 9. For {v,h) G Ct{M), define u^ by \2.10\) and let 7f be the associated 
flow map. Also define v"^ as u^ o rf . By 112. 12\) and 112. 13\) . we have 



sup ||Vo?r(t)||ffi.5(r) 

t£[0,T] 



,T||2 



^)) < C{M) 



(5.2) 



LL2(o,T;_H2.5(r) 

for some constant C'{M). 

We will make use of the following lemmas (proved in [TUl): 

Lemma 5.1. There exists Tq G (0,T) such that for all T G (0,ro) o,i^d for all 
V G Ct{M), the matrix a is well-defined (by 112.15]) ) with the estimate (independent 
ofve Ct{M)) 

\W\\L'^([),T-m(Q.)) + ||at||L~(0,T;ffi(a)) + Ikt II L2(o,T;H2(n)) 

+ ||aM||L=o(o,T;L2(o)) + ||ati||L2(o,T;Hi(o)) < C{M). (5.3) 
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Lemma 5.2. There exists Ti G (0,T) and a constant C (independent of M ) such 
that for all T e (0, Ti) and v G Ct{M), for all <j) e H^ifl) and t e [0, T] 

CMl^^n) < J^[\< + \DrM']dx (5.4) 

where 

In the remainder of the paper, we will assume that 

0<r<min{To,ri,f'} 
for some fixed T where the forcing F is defined on the time interval [OjT]. 

6. Preliminary results 

6.1. Pressure as a Lagrange multiplier. In the following discussion, we use 
H^'^{fl; r) to denote the space H^{n) n H^{T) with norm 

ll"llHi;2(f7:r) = II^IIhi(o) + Il"llff2(r) 
and Vy (Vc(T)) to denote the space 

Lemma 6.1. For all p e L'^(fl), t e [0,r], there exists a constant C > and 
4> e 7Ji'2(J7;r) such that al{t)(j)'j ^p and 

HWm-.^in-T) <C\\p\\L2(n)- (6.1) 

Proof. We solve the following problem on the time-dependent domain n{t): 

div(0o77(t)-i) =;po77(t)~^ in T]{t,n) := n(t). 

The solution to this problem can be written as the sum of the solutions to the 
following two problems 

diY{(f) o ri{ty^) = p o 'q{ty^ - p{t) in i]{t,n), (6.2) 

div(0o7y(t)~i) = p(i) in 77(t, f7), (6.3) 

where p{t) = — — ■ / p(t, x)dx. The existence of the solution to problem (|6.2p with 

l"l Jn 
zero boundary condition is standard (see, for example, [16j Chapter 3), and the 
solution to problem (|6.3p can be chosen as a linear function (linear in x) , for 
example, p{t)xi. The estimate (|6.ip follows from the estimates of the solutions to 
(lOl). D 



Define the linear functional on 7J-'^'^(ri; r) by (p, a^(t)(y5*)2^2(Q') where (^5 € H^''^ (yL\V) . 
By the Riesz representation theorem, there is a bounded linear operator Q{t) : 
L^[n) -^ iJi'2(r2; r) such that for aU ip £ H^-^{n; T), 

Letting ip = Q{t)p shows that 

\\Q{t)p\\m:2^n;r) < C\\p\\mn) 
for some constant C > 0. By Lemma \6A\ 

Iblli2(n) < \\Q{t)p\\H^-^n;r)y\\H^-^^n;r) < C\\Q{t)p\\Hi:2(n-r)\\p\\L''m 
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which shows that R{Q{t)) is closed in H^'''^{n;T). Since %{t) C R{Q(t))-^ and 
R{Q{t))^ C K(i), it follows that 

H'-'^in; r)(i) = R{Qit)) eHi,2(0;r) Mt). (6.4) 

We can now introduce our Lagrange multiplier 

Lemma 6.2. Let C{t) e H^'^{Q; T)' be such that C{t)ip = for any cp € V„(i)- Then 
there exist a unique q{t) G L'^{fl), which is termed the pressure function, satisfying 

\f^eH'-^\n;T), C{t){^) = (g(i),a^(t)^V)L^(o). 

Moreover, there is a C > (which does not depend on t (z [0, T] and e and on the 
choice of V £ Ct{M)) such that 

h{t)\\L-^(Q.) < C\\C{t)\\Hi;2(^-T)'- 

Proof. By the decomposition (|6.4p . for given a, let ip — vi + V2, where vi € Vy(t) 
and V2 G R{Q{t). It follows that 

for a unique ip{t) G R{Q{t)). 

^From the definition of Q{t) we then get the existence of a unique q{t) E L^{n) 

such that 

y^eH^-'\n;r), C{t){^) = (g(i),a^'(t)^:^.)L^(n). 
The estimate stated in the lemma is then a simple consequence of (|6.ip . D 

6.2. Estimates for a and /i. We make use of near-identity transformations. The 
following lemmas can be found in [9J and 10]. 

Lemma 6.3. There exists K > 0, Tq > such that if < t < Tq, then, for any 
{v,h)eCTo{M), 

||a^ - -^c'llL-(o,T;C»(ao)) - ^^ 5 (6-5^) 

\\a - /rf||L~(o,T;ff2(o)) < KVt ; (6.5b) 

||at - at(0)|U~(o,T;ffMo)) < C{M)t ; (6.5c) 

\\at\\L'^{o,T;m{n)) < K. (6.5d) 

We also need the following 
Lemma 6.4. For any {v,h) e Ctq{M), 

WHh^-HT) < CMt^'^ (6.6) 

for alio <t<To. 
Proof For {i},h) G Ct{M), \\h\\l,^^) + \\ht\\l2^r) < M. By h{0) = 0, 



\h{t)\\H2(T) < I \\ht\\m{r)ds < VMt. 
Jo 



/o 
Finally, the interpolation inequality 

Wylfimm-HD < C\\\7tf\\%\r)\\^lf\\ii\rr (6-7) 

implies 



Ii/3.5(r) < C\\hfj/!^^^\\h\\]i^^^.^ < CMt''\ 



D 
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Corollary 6.1. ||Li(i)||^i.5(p') and ||L2(i)||/fi.5(p) converge to zero as t — > 0, 
uniformly in {v,h) G CTo(Af). Furthermore, for t <1, 

\\L,mm-Hr) + WL^mm-^r) < C{M)t^/\ 

By the fact that ||/it||^2(r) < M and \\hu\\l2(^Q,T;H0-^^{r)) ^ ^^ if (^'M G CriM), 
similar computations lead to the following lemma. 



Lemma 6.5. For all (w, h) e CriM), 

||^t(i)||Hi.5(r)<CMii/« 
for alio <t<T. 



(6.8) 



7. The linearized problem 

Suppose that {v, h) S Ct[M) is given. Let f}{t) = Id+J^ v{s)ds and a = {\/f})~^. 
We are concerned with the following time-dependent linear problem, whose fixed- 
point V ^ V provides a solution to (|2.16p : 



vl-i^[a''eDi^{v)l],k^- 


-(5 


i 9),^ + F' 




in 


(0,T)xfi, 


(7 


.la) 


~a{v\j =0 in (0,T)xfi, 


(7.1b) 


\vDf,{v)\-qb'^a\Nt = a^ 


C~^{h){~\/oh,l) 


o?~f on (o,r) X r, 


(7.1c) 


-here 


\MCh){-^oh,l)]ofi^' 






ht o if = [h^^ o f)^]w„ - V, on (0, T) X r , 


(7.1d) 


V — uq on {t = 0} X f7 , 


(7.1e) 


/i = on {t = 0} X r . 


(7.1f) 


where Dfj{v){ = &fv^,^ + h^^v]^^, 6 = det(Vor)^), and 

1 


n 




with 


\/det(.go) 


y/detigo)A''^'''h^p 


,7<5 




^""'' = J/tV^^^^^'^) 


Gl'' _ (_!)«+- det(G,-J-i(l - 5„,)(1 - ,S^,)/i, J 


,a 




X (Gf - JK'Gi'Gth.K) 




and 



















Mih):^ ^det(G^)o7y^ Lf^{yXDh,D'^h)h,^0^+L2{y,h,Dh,D'^h) 

Here the thickness e is assumed to be 1. 

We will also use Lj^{h) to denote Cj^lh) -h M{h). 

Remark 10. £^ is a coercive fourth order operator for small h < 5. Actually, it is 
easy to see that L^ is coercive at time t = 0, and the coercivity of £^ for t > (but 
sufficiently small) follows from the continuity of h in time into the space H^ (F) . 
Moreover, by Lemma \6.4\ we have the following corollary. 

Corollary 7.1. There exists a ui > and < T <Tq such that for allO <t <T, 



^i||V^/(i)|| 



LHr 



< 



A''0'>'f,^l3it)fM{t)dS. 
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for all < t < T . Later on we will denote the right-hand side quantity of this 
inequality by Ej^{f), where the subscript h indicates that A is a function of h. 

Remark 11. Given {v,K) E V^{T) x T-t(T), for the corresponding ff , we have 

II^"'IIl=°(0,T;H2-5(j^)) + ll^t"llL2(0,T:ff2.5(r)) < C{M) 
where h2.13\) and i2.12\) are used to obtain this estimate. 



The solution of (j7.ip is found as a weak limit of a sequence of regularized problems. 
Definition 7.1. (Mollifiers on F) For e > 0, let 

RP := (1 - eAo)-i : H'{T) -^ H'+p{T) 

denote the usual self-adjoint Frederich mollifier on the compact manifold T, where 
Aq is the surface Laplacian defined on T . 

By the Sobolev extension theorem, there exist bounded extension operators 

Es:W{n)^H''{W), s>l. 

For fixed (but small) e and ti > 0, let p^ be a (positive) smooth mollifier on R". Set 
V = pe * Ei{v), F = p^ * E2{F), uq = p^* _E3(mo), where * denotes the convolution 
in space, and h = K^(h) for large enough m. Define fj and a in the same fashion 
as fj and d. Note that v ^ i e V{T), F ^ F in V'^{T), uq -^ Uq in H'^-^{n) and 
h->hin H{T) as e -> 0. 

The regularized problem takes the form 

vl-iy[a'^Df,iv)l],k = -iaU),k + F' in {0,T)xn, 

a\v\^ =0 in (0, T) X rj , 

^aJvq{-^Jiof]\\) + Kl\lv on (0,T)xF, 
htoff = \(h^^) o Tf\va. - V, on (0, T) X F , 



(7.2a) 
(7.2b) 

(7.2c) 
(7.2d) 



where 



^tif) 



V ^Uq 

h = 

e 



on {t = 0}xn, (7.2e) 
on {t = 0} xF, (7.2f) 



V^et{go) 



det(go)A"'5^V,a/ 



,7^. 



077 , 



Ml' = e 

h 



r ap-y 



L"''^i-,h, Dh, D%)h^ai3^ + L2(-, h, Dh) 



°v^iy,t)- 



Note that 41 (/) +M'j^=e [L-^if) 

7.1. Weak solutions. 

Definition 7.2. A vector v e Vv{T) with vt £ Vv{Ty for almost all t e (0, T) is a 
weak solution of |7.^| ) provided that 

h^iifi'^ofj-^) (7.3a) 



(i) {vt,v) + I I Alt' : Dr,^dx + a \_ A"^^'h[^p 



+ {if'- o f)--") dS + n / Ani; • Ao(^d5 = (F, ip) - a{Ml' , ^)i 



,7(5 



(ii) v{0,-) =uo 



(7.3b) 
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for almost all i G [0, T], where (•, •) denotes the duality product between Vy(t) and its 
dual Vv{ty , and h is given by the evolution equation ^7.2d\ ) and the initial condition 



h{y,t) 



hAv^ s)v"{rr^{y, s),0, s) + v%fj-^y, s), 0, s) 



ds 



(7.4) 



7.2. Penalized problems. Letting 6 > denote tlie penalized parameter, we 
define wg (with also e and ei dependence in mind) to be the "unique" solution 
of the problem (whose existence can be obtained via a modified Galerkin method 
which will be presented in the following sections): 



(i) (wet,^) 



2 ,^ D^we : D^^dx + <j j^ A'^^^'h^^ 



h^aiv" °V '^) 



+ iip''ofj ^) dS + K Aqw • AoipdS + {-div]j,al^y)L2(n) 
(ii) v{0, •) = uo 



(7.5a) 



(7.5b) 



where (•, •) denotes the pairing between H^{fl) and its dual, and h in (|7.5b ) satisfies 
(|7.4p with V replaced by we- 



7.3. Weak solutions for the penalized problem. The goal of this section is to 
establish the existence of v to the problem (|7.2p (or the weak formulation (|7.3p ). as 
well as the energy inequality satisfied by v and Vt- Before proceeding, we introduce 
variable qo and wi as follows: let go be the solution of the following Laplace equation 



Ago = Vuo : (Vuo)* - divF(O) 
go = i^{J^eiuo)lN,Nj ~ ^A^o'(O) + kA^^o • A^ 
and w\ be defined by 

wi = i^A-uo - Vgo + ^(0). 
By elliptic regularity. 



in 57 , 
onF, 



(7.6a) 
(7.6b) 

(7.7) 



~ l|2 



9011^-1(0) 



and hence 



<C 



I"o|Ih2(o) + \\Fm\i-{iX) + ll-^o^Wllffo-cr) + l|A^tio|||,o.5(r) 



< C{M) 



N^ol!/f-2(o) + ||wo||ff4.5(r) + 



I^(o)lli2(o) 



||^i|li.(f,) < C(M)[|l{toll'H2(,,) + |li*oll'H4.5(r) + WHmln^) + 1 

Remark 12. By i [6'.6]) . the constant C{M) in the estimates above can also be refined 
as a constant independent of M if T is chosen small enough. 



By introducing a (smooth) basis {ei)'^i of H^''^{il;r), and taking the approxi- 



mation at rank m>2 under the form wi{t,x) — ^ dk{t)ek{x) with 

fc=i 



hiiy,t) 



hAy^s)wf{r] '■(y,s),0,s) + ui|(?7 ^(y,s),0,s) 



ds, 



(7.8) 
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(i) (wm, f)L'-'(Q) + >^{alwu,j,a'^(p.k)L^n) + >^{{al°''i)tWi, V.k)L^n) 
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+ a / {A'^P'^%h\\,[~-h,{ip'^ o jy--) + ^- o ry-^^^.d^ 



(7.9a) 



+ a / A-'^-'^/i^^^^^ [-/.,,. (^- o ^--) + -h.v-{^% o f;--) + v-{^% o riryS 



^ {Ft,(p) - a ii '^ft.,a/37 + ^2 h,„{ip'^or] "■) - </?^ o 77 



dS* 



L?'3''/i 



Q/37 



L, 



htAv" ° r") - Kv^iv'^. ° ^-") - ^^"«« o ^-") 



dS' 



V (^ e span(ei, • • • ,e<?) , 
(ii) w«(0) = {wi)i, wi{Q) = (uo)£ in ^ , (7.9b) 

where qe = qo — -j^o^lw^j, and (uo)^ denote the respective iJ^'^(Jl;r) projections of 
uq on span(ei, 62, • • • , e^). 

Remark 13. The existence of Wk follows from the solution of 

d'l{t) + d'^{t)Aki{t)+di{t)Bki{t)+ I dfis)Cke{s,t)ds^F{t) 



for functions A, B , C and F ; however, the existence of the solution dk does not 
immediately follow from, the fundamental theorem of ODE due to the presence of 
the time-integral. A straightforward fix-point argument can be implemented, whose 
details we leave to interested reader. 



The use of the test function ip — wgt in this system of ODE gives us in turn the 
energy law 



1 d 
2Jt 



+ v{(d\a!l)tWi,^j,wu,k)L^(n)+v / (a^af)tw^ w^^ ^dx + K||AoWft| 



2 
L2(r) 



+ (g«,aiu;l^^.)L^(o) - {q,,-a'uw\t,,)LHn) - | I {A'^^"''Wl^phr,^^,dS 



-a {A^P-'%h\\., \hut + htA< ° V'l] '' dS + a A"P'<'hl\^x (7.10) 



-htA<t°v ") + V«"«t,«°^ l + v"iwkn°fi ") 



ei 



,7(5 



dS 



{Ft, Wit) - o- / {Lih,a(ii + -^2)(-Voft., 1) • {wit o ri '')dS 



-a/37. 



a {L'^"^h,^p^ + L2)v^ --hA<t,.°n ") + iwknOfj -) 



dS. 
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For the tenth term (the integral with — as its coefficient), we have 



By ei-regularization and the identity 



/r Vdet(5o) 



^f^:^[^)[A^^^% 



,7(5 






4A^{^){A^f"^\ 



IT y^dcUso) 






we find that 



{A-^'")thl^^h\l,^^,dS 



< C(ei) 1 + ||/it||//2.5(r) ||Vo/i£||L2(r) \\wi\\m(n) + \\wu\\m(n) 

Similarly, the second part of the eleventh term and the last term of the left-hand 
side can be bounded by 

C{ei)\\ht\\m-Hr)\\^lht\\mr)\\wu\\m(n) 

where we also use the ei-regularization to control V^wu- It also follows that the 
last two terms on the right-hand side can be bounded by 

C{M) l + \\ht\\H^.^T) \\wu\\m(n)- 

With positive 9, the fourth term of the left-hand side involving the square of qu 
acts as a viscous energy term. Integrating ()7.10p in time from to i, we then get 

ll^^ttlli^ro) + l|Vo/ift|li2(r) + / \\^wu\\l2u^)+ K\\wtt\\]i2rp)+ S\\qtt\\l2n-i) ds 

Jo ^ J 

< C{M) \^\wam\l^n) + 1 1 "'^ (0)1 1 HI (o) + lk^(0)ll?,o.5(a)] (7.11) 

+ C(ei) j \l + \\Ms)\\m-.iT)]\\'^lhH{s)\\l2(r)ds 
Jo L J 

Jo Jo ^ -' 

where C(ei),C(6') — > oo as £1,6* ^ 0, and we use 

ll/(i)l|x < ||/(0)||x + f \\It{s)\\xds < ||/(0)|U + Vt f WMsWxds 
Jo Jo 

for / ~ Wi, f = he and f — gg to obtain (|7.1ip . 

Remark 14. The 9 -dependence follows from estimating the terms {qu, dlfW\ j)L'^(m- 

1 
29' 



(gft,a^jW^j)i2(o) 



..i Il2 



< -l|'?ftllL2(o) + :j^\\4t\\l^{n)\\Wij\\l2(^n) 



< 



C{M) 



llfelli2(o) + ^^[||V«;,(0)||i2(o)+iy l|Vu.,t||i2(,,)(s)ds 
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Jo '- 
+ 0\\qit\\l2^n)]ds<Cie,,9)Noiuo,F) 



(7.12) 



where 



No{uo,F) := ||tio||H2.5(o) + l|Mo||ffi-5(r) + ll^t|li2(o,T;ffi(n)') + \\P{0)\\Ho^n) + 1- 

We can then infer that we defined on [0, T], and that there is a subsequence, still 
denoted with the subscript £, satisfying 



where 



wu^wgt inL2(o,T;iJi'2(r!;r)) 
Vlhi^Vlhe mL^iO,T;L^{T)) 



qe=qo~ n^-iwlj. 



(7.13a) 
(7.13b) 
(7.13c) 
(7.13d) 
(7.13e) 



From the standard procedure for weak solutions, we can now infer from these weak 
convergences and the definition of wg that wut G L,^{0,T;H^{ny). In turn, wu E 
C°{[0,T];H^{ny), Wi e C°{[0,T];L^{n)) with wg{0) = uo, wet{0) = wi. 
Moreover, (|7.13p implies that wg satisfies 



(i) / iwgtt,v)L^n) + v{alwgtj,a^^^k)L^{n) + v{{a{a^)twg, (p,k)L2(n) 



dt 



aia'^wlt,ip\dx + v I {ala'l)twl.ip'^kdx 



Jr 



(7.14a) 



T 



+ K / AoWgt-AoipdSdt- / {{a'iqg)t,(pj)L^(n)dt 



"T 
T 



UFt,ip)-(7 L"^'^/l^a/37 + ^2 /l^cr ((/?"■ 077 ■^) - y)^ 077 



dS' 



L^^'^T/i., 



■,a/37 + -Zj2 

-v^{ip%ofj-^)Yds} 
(ii) wet(O) = w)i,we(0) = uo in f7, 



/l,,,(V.-ojy--)-/l,,i;-(^-„o^--) 



(7.14b) 
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for all (f e L?{{),T\H^''^{Vt\V)). Choosing (p to be independent of time, wc find 
that for aUte [0,r], 

iwet,f)L^(n) + o / Dfiiwe) ■ D,-^{ip)dx + k Aowe • AoipdS 
^ Jn Jr 



{F,^)+<jJ [Lf'<'h^^f,, + L2 



h^cr(p'^ OTj "^ + (p^' orj 



dS + c{(p) 



for all ip e H^''^{Q; T), where c{tp) G M is given by 



c{ip) = {wi,(p)L^(m + o / Def(uo) : Dciifdx- (go - - divuo, div(p)L2(o) 
- (#(0), (p)i2(n) - aiM'o' (0)(0, l),ip)L2ir) + k(AoUo, AoV5)L2(r) . 

By compatibility conditions (|7.6p and (|7.7p . c((y9) = 0. Therefore, the weak limit 
(we, he) satisfies, for all t e [0,T], 



(wet,</')L2(n) + 77 / Df^{we) : Df^{ip)dx + k AqWo ■ AoipdS 



{criqe,^,)L^(n) + o \ A"^^^/j^;„^[-/i,.(^- o ff^) + ^^ o ^-^];^,d5 (7.15) 



{F,^)-G 



^•1 n^ap-f + J^2 



h^^ip" ori '^ + (p^ OT] "■ 



dS, 



iovallipe ffi'2(J7;r). 

Since we € L^(0,T; _ff^'^(r2; F)), we can use it as a test function in (|7.15p and 
obtain (after time integration) 



2ll^0lli2(o) +2^hih'0'] 



2 ll^l^elli^Co) + «ll Aowe||i2(r) 



+ ^lke|li2(j^) ds-9 I {qg,qQ)dt 



2 7o 



(^"''"')*/^«!«0^^:..'^^'^^ (7.16) 



Consequently, 

lhe(t)||i2(a) + ||V2/i^^(i)lli2(r)l + / WVweWh^^^ds + n f \\we\\l2(^r)ds 

^ Jo Jo 

k0lli2(o)rfs 

<C(M)[||uo||i2(o)+0||9o|li2(o) + ll^lll/l(o)' + ll-^|^-Voft°ry^l)||i2(^) 

+ C(Af) / ||/i,||^2.5(r)||V2/i^^||2 ds 
Jo 
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where 

Ni{uo,F) = \\uo\\H2(n) + ll"o||ff4.5(r) + \\F\\l2(^Q^T-HHny} + ll^tllL2(o,T;ifi(iy') 



WFrnwinn) 



By the Gronwall inequahty, 



sup 

0<t<T 



\we{t)\\l2(^a) 



rT , 



-■ ./n L 



<CiM)N,iuo,F). 

7.4. Improved pressure estimates 

as, for a. a. t e [0,r], 



ds 
(7.17) 
. By ei-regularization, we can rewrite (|7.15p 



^ Jn 
+ a I C)l {hl')[ - h^ o fj-'ip- + <^-] dS = (F, ^) + a(>l^^i (-Vq/i o f]\ 1), (p)r. 

Therefore, by the Lagrange multiplier lemma, we conclude that 

Ike|li2(!:2) < C{M) WwetWHi^ny + W'^weWh^n^ + \\F\\jji^^^y + k|| Ai^w0||^-2(r) 

and hence 

lke|li2(o) < C{M) [WwetWh^n) + W^wgWl^^,,^ + K\\we\\jf2^r) + llVg/ie|li2(r) 

+ \\F\\Jii<^ny + '^]- (7.18) 

7.5. Weak limits as 6* — > 0. Since wgt & L^{0,T;H^''^{Q;T)), we can use it as a 
test function in ()7.14p . Similar to the way we obtain (17. lip , we find that 

^WwetWh^n) + 2 / \F>v^et\\l2i^^jds + -Ef^{hl\) + k WAlwetWhf^r) 

+ ^ I ket 1 11,2(0)^5+ / iqet,al^wl.j)L2(n)ds- / iqg,dlwlt^)ds 
Jo Jo Jo 

■,C{M)No{uo,F) + C{M) f Mt')\\hin) I \\^wet{s)\\l2^n)dsdt' 

Jo Jo 



< 



ds 



+ C(ei) / l + \\ht\\H2.^r) W^lKWll^^^^ds 



Bydng) 



< 



iq0,alwlt,j)ds <C{M,S) heWh m)ds + S \\Vwgt\\l2m)ds 

Jo Jo 

C{M) [7Vi(wo, F) + j (\\wet\\l2(n) + i^WweWmiv) + l|Vg/ie||i2(r))ds 



+ 5 liVwet||i2(o)rfs 



(7.19) 



where (|7.17p is used to bound ||Vwe||^2(o T-L^fn))- 
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Integrating by parts, 



t 
{qet,aJtwlj)L2(n)ds = (ge, a-tW^j)L2(n)(i) + (go, <»"oj)i^(n) 



;-,ds. 



By e-reg 



- / (g0,a^ttW^j)L2(n)rfs- / (g0,a^fU;^tj)L2(o)C 
Jo Jo 

;ularization, the last two term can be bounded by 

C(e)||Vwe||L2(o) + \\^W0t\\L^n) ds 



C{M) f ||ge||L2(n) 



and hence 



) Jo 

,rfs + (5 / ||Vwet|li2(o)ds 
Jo 



For 



<C{M,S) f \\qe\\Un)ds + C{e) f \\S/we\\h^n) 

Jo Jo Jo 

<Cie,6)N,iuo,F) + CiM,d) f WwetWha^.ds + C'iei) f 

Jo Jo 

+ 5 \\Wwet\\l2m)ds. 
Jo 

{qe,dj^Wg,j)L2(n){t), it is easy to see that 

2 ' ' '' ' 



\\^lhe\\l2^r)ds 



(7.20) 



^hile for {qo,i 
(fr20| and the 



while for {qo , Uq ^Uq j) l^ (n) , it is bounded by C{M)Niyuo, r j. v_/umuimiig 
/Pt-ot^ „„j 4-u., estimates above, by choosing S > and (5i > small enough. 



\\wet\\h^n) + \\F\\LHn) + l^ 
\{uo,F). Combining (1719)) 

1 Ai ^ n Qmnll pnnncrVi 



" ' l|Vo^0tllL2(r) 



/■* r 



""^6*11^2(0) 

iV2(uo,J^)+ / {\\wgt\ 



<C(ei,e) 



l|2 

IIl2(s-7) 



l|Vt«et|li2(n) + K||w0t||^2(r) + 6'lket|li2(o) 
l + ||/it|U2.5(r))||V2/ie,||2^(P) 



ds 



fIIh3(o) 



where N2{uo, F) = iVi(uo, F) + ||-F||^oo(o T-L^(n))- ^y ^^'^ Gronwall inequality, 



ll^etlli2(j^) + ||V^/iet|li2(r) 

/o 



|Vwet||^2(o) +K||wet||^2(r) 



ds 



Jo '- 
< C{ei,e)N2{uo,F) + Ci{ei,e)\\^wg\\l2^n)- 

By using wg{t) — uq + wgtds, we see that 
Jo 

!|Vwet||i2(o) +K||wet||H2(r) ds 



(7.21) 



lk0t|lL2(n) 



||Voft,et||^2(r) 



< 



C(ei,e)7V2(uo,F) + Ci(ei,e)i / \\Vwgt\\l2,n.ds. 

Jo 
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Therefore, for any < t < ti = niin 



H--^l 



we have 



1 ft r 



ds 



<Ciei,e)N2{uo,F). 
By wg{ti) = mq + / wgtds, we also have 

II Vz«,(ti)|ii.(^) < C(ei, e)7V2K, F). (7.22) 

For t > ti, shice we{t) — wg{ti) + / wgtds, we have from (|7.2ip and (|7.22p that 

WwetWht^n) + l|Vo/iet||i2(r) + / llVweiHisjo) + ^Ihetl' 



l_H"2(r) 



ds 



<C{e,,e)N2{uo,F) + Ciie,,e) \\wg{h)\\l,^,,) + {t - h) / ||Voi«et||i2(f,)ds 



<C(ei,e)7V2(Mo,J^) + Ci(ei,e)(t-ti) / \\'^oWet\\i2^a)ds 

Jti 



Therefore, for any ti < t < 2ti, we also have 



1 



2 7o 



Wweth^n) + l|Vo^9t|lL2(r) + 9 / l|Vwei||L2(o) + K||wet||i/2(r) 



ds 



<C{ei,e)N2{uo,F) 



2ti 



which with wg(2ti) = ua + wgtds gives 

Jo 

||Vu;e(2ti)||i.(,,)<C(ei,e)iV2K,i^). 

By induction, for any t G [0, T], 

\\wet\\h(n) + l|Vo^et||i2(r) + 2 / [ll'^«^et|li2(o) + nWwetW^i^^r)]^^ 

<Cieue)N2{uo,F). (7.23) 

We also get a ^-independent bound for ||g6i|||2(o T-L^(n)) ^^ (|7.18p : 

lk«lli^(o,T;L^(o)) < C(ei, e)iV2(wo, i^)- (7.24) 

Let 6* = ^. Energy inequalities (|7l7l) . (rr23| and (pOil show that there exists 
a subsequence w j_ such that 

mg 

w 1 ^ D 



in 2.2(0, T;i?i'2(f7;r)) 

w^t^Ot in L^{0,T;H^''^{n;T)) 

Vo/i_L^Vot) in 2.2(0, T;i2(f])) 

mg 

Vlh^t^Vll)t in l2(o ^..^2^^)) 

<?_L^q in L2(^0,T;L2(f])). 



(7.25a) 
(7.25b) 
(7.25c) 
(7.25d) 
(7.25e) 
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Moreover, (|7.17p also shows that ||a^K;*j_ -\\L'^{o,T;L^(n)) — > as m ^ cx). Therefore 
the weak hmit satisfies the "divergence-free" condition (|7.2l b'). i.e., 

e Vs(T). (7.26) 

Since (|7.17|) is independent of 6 and ei, by the property of lower-semicontinuity of 
norms, 



sup 

0<t<T 



P^"llL2(0,T;L2(n)) + '*l|0||//2(-p) 



\m\\hin) + \\^lm\\h(r) 

<C{M)Ni{uo,F). (7.27) 

By (|7.25|) and ei-regularization, the weak hmit (o,f),q) satisfies, for all (p E 
L^iO,T;H^'Hn;T)), 

/ {Vt,ip)L^n)dt+^ j Df,{\:i) : Dfj{ip)dxdt + K Aq\:i ■ AoipdSdt 

Jo ^ Jo Jn Jo Jr 



{alc\,(p]j)L^(n)dt + cr 



T 

Jr 



dS \dt . 



By the density argument, we find that for a. a. i G [0, T], </? e H^''^{Vt\ F), 



(Ot,(p)L2(o) + TT / D^{v) : Dr,{ip)dx + n I AqO • Ao(pdS' - (a^q, (p^)L2(n) 



+ a / A^P^%^l^J-h,,{^'^ o f,-^) + ^- o rl^'sdS 



(7.28) 



(F, ^)-a 



- afijS 






h_cLp'^ o T] '^ + ip^ o rj "^ 



d5. 



or after a change of variable y' — t]'^ {y, t), 



{Vt,(p)L^{n) + 7:{Df^t>,Dfjip)L2(^^^ + K AoD • Ao<^d5 - (a:^q, vj^Oi^cn) (7.29) 



+ a C}At))i-Vohofj^,l)-ipdS = {F,ip)-a M}^ {-V oh o fj\ 1) ■ ipdS . 
Furthermore, if (p G V5, then 



(Ot, 'yj)L2(o) + ^7(^70, Dfjip)L2i^^^ +K AoO • AoipdS 



+ aj C';^m--Vohorf,l)-ipdS^{F,ip)-aJ M'^i-^oho f ,1) ■ ^'^dS 
for a. a. i G [0, T]. In other words, (o, (), q) is a weak solution of (|7.2p . 

8. Estimates independent of ei 

8.1. Partition of unity. Since Q. is compact, by partition of unity, wc can choose 
two non-negative smooth functions ^0 and ^1 so that 

Co + Ci = 1 in f^ ; 
supp(Co) CC n ; 
supp(Ci) CC F X (-e,e) := ^i. 
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We will assume that Ci = 1 inside the region n[ C fli and Co = 1 inside the region 
n' C n. Note that then Ci = 1 while Co = on T. 

8.2. Higher regularity. 

8.2.1. ei -independent bounds for q. Similar to (|7.18p . we have 



\F\ 



L^(n) 



(8.1) 



8.2.2. Interior regularity. Converting the fluid equation (17. 2p into Eulerian vari- 
ables by composing with fj~^, we obtain a Stokes problem in the domain fj{n): 



-i/Au + Vp = F o ?7 1 - Ot o ?7 1 + i/a^_j. o -q ^u,^ - pa^^^ o V ^ , (8-2a) 

divu = 0, (8.2b) 

where u ~ vofj^^ and p = qofj^^. By the regularity results for the Stokes problem, 



< C 



\\P°v TL^f-m) + II"* °^ ^WUvm + II^HIi^wfi)) + llPlli^(f;(o)) 



^ll_H"i-5(r) 



or 



I|o||ff2(o) + Ilq||?/i(i2) < C[\\F\\l2^^) + lloJi^fo) + l|o|lHi-5(r) 
for some constant C independent of M, e. By (|8.ip . 

I|o|Ih^(o) + hWrnm < GW [llfJi^fo) + l|vo||i.(f2) + iion^^.^p 



72(,eil|2 



+ l|i^ll 



+ ll^o')''llL2(r) -r iKiiL2(jj) 
Similarly, 

ll'^ll_f/3(o) + ||q||_H-2(o) < C* ||-F||/i-i(o) + l|0t|li/i(s:2) + l|o||_H-2.5(r) 

+ C{M)[\\\/X>\\l^^,,^ + \\^\\j,^^^^ 

and therefore by ((57T|l and (jO)) . 

l|o||?f3(n) + Ilq||ff2(j,) < C(A/)[||tiJ^,(^) + ||Vo||i.(f,) + llV^oH^^j^^) 

+ l|Vgf]^MIi2(r) + 11^111/1^7) + 1 



(8.3) 



(8.4) 



For the regularized problem, because the e-regularization ensures that the forcing 
and the initial data are smooth, while the ei-regularization ensures that the right- 
hand side of (|7.2c[) is smooth, by standard difference quotient technique, it is also 
easy to see that 



VaV e L^{0,T; H\ni) n H^{r)) for fc = 1,2,3,4 



(8.5) 
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7) 



in Q 



Since (|7.25b ) implies that 0* E L'^{Q,T;H^{il)), by ei-regularization and ((Ki)) we 
conclude that 

e L2(0, T; H^in)), q e L^^q^ T; H^{n)). (8.6) 

8.3. Estimates for Of(0) and q(0). By (18. 6|) and ei-regularization, (t),f),q) satis- 
fies the strong from (|7.2p . Taking the "divergence" of (|7.2a[l and then making use 
of condition (|7.2l b). we find that 

-^O^fc - mf [ap^(t))^] jfe = -af (a^q),,fe + a.^F^ 

Let t = 0, by the identity a^^ = —al.v\a^j, 

Aq(0) = Vuo : (Vuq)^ - div(F(0)) 
with 

q(0) = i^{Dciuo)iN,Nj - aMo'{0) + kA^uq 
while ()7.2a|) gives us 

Ot(0) = i/Amo- Vq(0) + ^(0) in Q. 
By standard elliptic regularity result, 

l|f*(0)|li.(o) + l|q(0)||^i(o) < CNo{uo,F) (8.8) 

for some constant independent of M, e and ei. 

8.4. LjL^-estimates for Ot- Since Vt G L^{0,T; H^{i^)), we can use it as a test 
function in (|7?29l) . By (fr26l) . we find that 

llOtlli^Cn) + ^^ / lAiOpda: - ^ j (7^^o)^a^^tO'fc'^a; + k J AqX, ■ Ao^dS 



on r 



q4tO}da; + cr / C'jl{\)){-ynhoff,l)-X)tdS 



{F,x>t)-a M'j^{~Vohofj^,l)-tttdS. 



By (ESI, 



/ iD^v)ld^,^,dx < CiM)CiS)\\Vx,\\l. 
Jn 



(O) 



s\y^ 



H^{n) 



and by (|8.ip and the interpolation inequality, 



e ^k 



qa^o,^ 



^dx 



< CiM)C{S) [||Vo||i.(o) + llV^f)^^lli.(r) + |lF|li.(,,) + 1 

for some C{6). Also, the last term on the left hand side is bounded by 

CiM)[\\Vti)''\\LHr) + l]\MHHn) 

<CiM)C{S^)[\\Vti)^^l.(r) + l\+Si\MJfi<^ny 
Combining all the estimates above, 



2ll"*llL2(o) 



< C 



^ d [ ,^ ,P , K. d 

/ \Df,<ordx^ 

4 di it^ ' ** ' 2 rfi jr 

74[,ei||2 



A.vVdS 



W^Hi^n) + llV^f)^^ lli.(r) + ll^lli2(o) + 1 + 5\\x>rH.^n) + '^illodll/Ma) 
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for some constant C depending on M, S and 6i. Therefore by (|7.27p . 



/ I|0t||i2(fj)ds + ||VD(t)||i2(j^) + K||o||^2(r) 
Jo 



< C 



N2{uo,F) + 



\<h''\\h(r)ds 



sj Ml2(^n^ds + 6iJ 110 



8.5. Energy estimates for VqI; near the boundary. Because of 
in (|7.28p can be used as a test function in (|7.29p . It follows that 
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(8.9) 
t\lH^(n)d-s- 
,Vg(C?V2t,) 



■^n 



<C{M,d3 
By (U^D, we find that 



<C{M)\\\Vlt)''\\H^r) + l 
1 



"ll-ff^ir) 



2 

LHi(r) 



'^3||o||_H-4(r) 



||f)||^4(r) <C(e) / ||/i||^5(r)||D||//4(r)ds < C(e) / ||o|||^4(r)ds 



and hence 



el/'hel^ 



4Hf) 



< c 



■^K 



LH"*(r) 



(-Vofto^^,l). V^od^ 



<53l|0"" 



_ff4(r) 



for some constant C depending on M , e and 53. Since 



Ao/ 



d 



y/d^ti^gfg^f 



by the regularity on T (and hence on go), 

J lAoV^ol^d^ < y A^D . {Vt>v)dS + CMhhdMhht) 

<J^^lv i^tv)dS + C(<5)||t,||^,(,,) + <5||t,||^.(p) 
which implies, by choosing (5 > small enough, that 

^2|M!ff4(r) < y A2« . {W^,v)dS + C\\v\\l^^,,y 
By the identity 

= (q, V^ai (C? V^o'^),,) +4(CiVoq,Vo4Ci/V2o'=) + 2(Voq,Ci'VoaI.V^t.',) 
- 2(CiVoq, VoiaiCi./Vy) + 2(q, VoiaiCi.i^oCi^y)) (8.10) 

+ (Voq,Vo(CiVo4Vot';^£)), 
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(|0)) and ([13]) imply that 

(q,ai(V^'(CiVgo'^),£) < C{M)\\q\\H^n)\MHHn) 
< C{M)C{S) \\Mh^n) + l|Vt.||i.(,,) + II VVoO||i.(,,,) + kMJ,. 



(r) 



'''"lL2(r) + II-^IIl2(s;2) 



SMls 



(n)- 



For the viscosity term, 
Jn 

J Si 

Vo(a,^a-)VoDf^ + Vo(af a ■)VoO;^] (Ci Vgt)^),fcda; 
and hence by interpolation 



C{M)C{5) ||Vti||i.(j,) + ||VVoti||i.(o;) + 6\Ml^ 



m- 



Summing all the estimates, by letting ^3 



ld_ 

2di 

< C 



V , 



ClVoO||i2(s-j) -r ^IM^Jl voiJ||£2(s-j) 



we conclude that 



-2«,|,|,2 



^tWht^n) + IHmm + l|VVot)||i2(o;) + l|o|lH2(r) + W'^'t^^'WUr) 



2 n-nH4(r) 
2 I IIV7V7 i,l|2 I ||t,||2 I llv72Lei||2 



1^1 



H^n) 



1 



C / ||o||H4(r)'^s + 5||D||^3(o) 
Jo 



for some constant C depending on M, n, e and (5. Integrating the inequality above 
in time from to t, by (|7.27p we find that 



ivguwr 



L2(ni 



|VV2o||i2(,,^) 



K f74 



//4(r) 



ds 



<CN2iuo,F) + C I ||oJi.(o) + l|VVoti||i.(f,,) 



+ C 



lH2(r) 



(8.11) 



Mr)\\H4{r)dr + 5 / ||o||^3(n)rfs- 



By using Vo(CiVot)) as a testing function in (|7.29|) . similar computations leads to 



l|Voo(i)llL(Oi 



(Oi) 



l|VVoo||i2(Oi)+K||o|lH3(r) 



ds 



< C{M)N2iuo,F) + C{M,S) 



t ^ 



mlmn) +'«l|0|l_ff2(r) 



ds 



^.12) 



Jo Jo Jo 
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8.6. Energy estimates for vt - LjiJ^-estimates. In this section, we time differ- 
entiate (|7.29p and then use Ot as a test function to obtain 



(titt,Ot) +v 



^eiD^^)lk 



Hldx + a 



4nr)(-Vo/io^^l) 



■X>tdS 



+ K / |AoOtrd5- / {aU)tX>1idx^{FuX>t)~cT / M'^{-^oho rf ,1) 
JT Jn ' Jr '- 

By the chain rule, 






■ VtdS. 



e 



[L^(r)]^HVo/i,-l)] orf-t^tdS. 



By using H^{T)-H ^(F) duahty pairing with e-regularization on O and v, it foUows 
that 



{C^{\)^^) + M'j^){-^^hof,^,\) 



■VtdS 



< C{e) [\\Vlt)\\L^r) + WylhtWmr) + ij IMhht) 
<C(e,d3) 



|f||ff4(r)ds + ||o||^2(r) + l 



S3\Ml2^r) 



< C 



|o||^4(r)rfs+l!ti||J^i(j^) + l 



SMH^in)+S3\MH2^r) 



for some constant C depending on Af, e, S and 63, where we estimate ||t)||^2(r) by 
interpolation. 

Also by interpolation 



/ \DfjVt\''dx^2 



< 2 



«'^^(fH 



0^ i^dx ~ 2 



{d\-a%Xi\ + {d!l-a])t'o\ 



^ikdx 



d'yDrM\^^ikd^ + CiM)C{S,S^)\\Vv\\l.f^^^ 

Note that 

{Ft,t,t) < CWFtWHiinylMmm < CiS,)\\Ft\\l.^^y + <5i||ot|lHi(n)- 
Summing all the estimates above. 



1 d 
2dt' 



^tWhi^n) + Tl|V0t||i2(o) + nWAoVtWh^r) 



< C 



||o||^4(r)ds + ||o||^i(j,) + 1 + C{6i)\\Ft\\j,.^^y 



(8.13) 



-S\Ml^n)+Si\MHnn)+h\MJi2(^r)+ / idU)tvledx 
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for some constant C depending on M, k, S and 6i. As in [TU] and [TT], the integral 
involving the pressure q has the following estimate: 



/ iaic\)tvl^dxds < C{M)C{5, 5i)Nz{uq, F) + 5 f 
Jo, ' Jo 

t 



0|lff3(n)rfs 



+ •^1 / Pt\\m(n)ds 
Jo 



where 



N3{uo,F) := |luollff2.5(o) + llMo|llf4.5(r) + ll^lli2(o,T;Hi(o)) 
+ \\Ft\\L^{o.T;m{ny) + I1^(0)II_h-i(o) + 1- 

Integrating (|8.13p in time from to t and choosing 61,63 > small enough, 
(rr27| and (gJl) imply that, for all t£ [0, T], 



iio*wr 



L^(n) 



|VOt||i2(fj)+K|!ot|i^2(r) 



ds 



<CN3iuo,F)+C / ||0(r)||^4(r)drds + (5 / \Mmmds 



(8.14) 



for some constant C depending on M, k, 6 and 62- In (|8.14p . (|8.8p is used to bound 



\\vmr 



.7. ei-independent estimates. Integrating (|8.3p in time from to t, (|7.27p . 
and ((51^ imply that 
t 



1^ I ii^ii^ 

Lf/2(n) + llMll_ffi(o) 



ds 



< CiM)N,iuo,F) + / ||t,,|li.(o) + Mhhd 
Jo '- J 

2 j„j„ , X I ll>,l|2 



ds 



<CN3iuo,F)+C / ||o(r)||^4(ndrds + (5 / ||o||^3(n)ds (8.15) 

Jo Jo Jo 

for some constant C depending on M , n and 6. Integrating (|8.4p in time from to 
t, making use of (|8.1ip . (|8.12p . (|8.14p . (|8.15p . and then choosing ^ > small enough 
and T even smaller, we find that 

^MJfs^n) + Mmin)]ds<CN3{uo,F)+C f f ||o(r)||l,.(p)drds (8.16) 
^ Jo Jo 

for some constant C depending on A/, k and e. 

Having (|8.16p . by choosing 62 > small enough, the estimates (|8.1ip can be 
rewritten as 

ds 



l|V^o(i)||i.(,j^)H 
<CN3{uo,F)+C 





t ps 



\Hr)\\ 



|^4(■p^uTuS 



5.17) 



for some constant C depending on Af , k, and e. Therefore, 

X{t)<c\ f X{s)ds + N3{uo,F) 
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where 



Xit) 



JQ 



By the Gronwall inequahty, 

nt nS 



"'0 



Hr)\\'\j,^Y)drds < CN:i{ua,F) 



(8.18) 



for aU t E [0,T] for some constant C depending on M, k, and e. Having (|8.18p . es- 
timates (|8.9p , (|8.14p , (|8.16p and (|8.17p along with the standard embedding theorem 
lead to 



sup 

0<t<T 



\^i^)\\H^(n) + ll>^t(0llL2(O) + ll'^llv3(T) + lllllL2(o,T;H2(n)) 

+ 4^\\hio^T;HHr)) < Cmiuo,F) (8.19) 

for some constant C depending on M , n and e. 

8.8. Weak limits as ei -^ 0. Since the estimate (|8.19p is independent of ei, the 
weak limit as ei -^ of the sequence (o, f), q) exists. We will denote the weak limit 
of (o, [), q) by (wk, h^, Qk)- By lower semi-continuity, (|8.8p and thus (|8.19p hold for 
the weak limit (v^, ^k, 9k). Furthermore, 



^ Jn Jr 



+ K AqWk • AoipdS ~ (^K, afc</'/)-L2(0) 
= {F,^)-aJ e[[M{h){-'^oh, 1)] o f]^ 
for aU (y9 e H^'^iQ; T) and a.a. i G [0, T]. 



[£s(/i„)(-Vo/i,l)]of/- 



lysdS' 



ipdS 



(8.20) 



9. Estimates independent of k and e 

9.1. Energy estimates which are independent of k. Although (|8.19p doesn't 
imply that /i^ G i7^(r), h,^ is indeed in _ff'*(r) by (|7.4p . Therefore, we have that 
(wkj/ikj^k) satisfies 

u«,j - i/[a^i:>^(w«)^],fc = -(afq^),fe + F^ 



a'v^\ = 



in(0,T)xrj, (9.1a) 

,^j „ in(0,T)xf], (9.1b) 

[^//^^(i;,)^' -g«,5^']aj^iV, = ae[£^(/i«)(-Vo/i,l)]ofy" on (0,T) x T, (9.1c) 

-I- aQ[M-h{~Voh, 1)] o 77^ + kA^i;^ 

/ito,7- = p,„)o^-]i;„-«, on(0,T)xr, (9.1d) 

V = uo on {f = 0} X ft, (9.1e) 

/i = on{i = 0}xr. (9. If) 

Having (pix)) . (|A.7p in Appendix El implies that /i^ is in i/^(r) for a.a. t G [0, T] 
with estimate 







t rt 



|Vo^K||i2(r) + lk«||^3(o) + lkK||ff2(n) + 



ds. 
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where the forcing / in (jA.7[) is given by 

[lyDrjMi - QJlHrn ~ ae[MfX-^oh, 1)] o f]\ 

By the same argument, (|7.18p holds with all 9 replaced by k. Therefore, by (|8.4p 
(which follows from (|7.18p ). 



t ^ 



\\^ohn\\imr)ds < C{e) / Wv^tWimn) + l|V^/i«||i2(r) + ||VXIIhi(Oi) 
Jo '- 

+ C{e)N2{u„,F). 



ds 
(9.2) 



With this extra regularity of /ik, the energy estimate (|8.19p can be made inde- 
pendent of K. In Appendix IB] 2. we prove that 



yl|V^/..(t)||i.(r) < 



e 



[L-^ih,)i-VoK 1)] o V" ■ yt,{Ct^t>v,)dSds 



+ C' 



"T 

r,l|2 , iir 112 , ||r,||2 



1 + Mmin) + I|/it|lff2.5(r) + \\h\\H^r) \No!^'^\\L^r)ds 



C 



'*llff5(r) + l ds + S \\vK\\H3n2)ds + 5i / \\Vlh^ 



{nr 



H3(r) 



ds 



for some constant C depending on M, e, 6 and Si. By (|9.2p . 



yl|V^/iK(t)||i.(r) < 



e 



Jr 



[ih(/^«)(-Vo/i, 1)] o 77^ • Wl{CfWlv^)dSds 



+ C'N2{uo,F) + C' \\VW\\l.^n,)+K{s)\\Vth.\\h^r) ds (9.3) 



JO Jo 



where 



K{s) := 1 + Mls^,,^ + WhtWU^r) + mlHry 
With (US]), (|5TT|) now is replaced by 



\\^lMt)\\hin,) + \\^th.{mh^r) 



t ^ 



\^^l'"d\^n^) + l^h^^H'^iT) 



ds 



t ^ 



<C'N2{uo,F) + C' / \Kt\\Un) + \\^W\\HHn,)+K{s)\\Wth^\\l,^r) 



ds 



+ S I \\vn\\m{n)ds + Si / \\vnt\\m{<:i)ds 



(9.4) 



for some C" depending on M, e, 6 and 6i, where (jA.Sp is applied to bound '«||wK||^3(r) 
(this is where ||wKt||^2ff2') comes from). Similar computations leads to 



Vo««(i)lli2(Oi) + llVil/iK(<)lli2(r) 

< CN2iuo, F)+C f \\Vthn\\l2(^r)ds + 6 f 
Jo Jo 

for some constant C depending on AI and S. 



lVVoWK|li2(ni) +K||wK||^3(r) 



\'"K\\H3(n)ds 



ds 
(9.5) 
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In Appendix [Cl we establish the following k- and e-independent inequality for 
the time-differentiated problem: 



ft rt 

Jo Jo Jr 

rt 



[Lj,{K){Voh,-l)]ofj^ 



VntdS 



+ CN3{uo,F) + C K{s) \\^th.\\i2^r) + mh.t\\i2^r) 



ds 



+ {5 + Ct^'^) \\v.\\l.(n)ds + iSi+Ct'^') \\v.t\\'HUn)ds + S2\\yth.\ 



2 



for some constant C depending on M, 6, Si and ^2- Therefore, (|8.14p can be 
replaced by the following estimate: 



ll«KtllL2(o) + l|Vo^Kt||L2(r) +/ l|VwKtllL2(o) +K||Aov«:tlL2(r) 





4i, l|2 



ds 



<CN3{uo,F) + C Kis) W'^th.Wh^D + W'^oh.trmr) 



ds 



(9.6) 



Jo Jo 

9.2. K- independent estimates. Just as in Section [5771 we find that 

<C{M)N2{uo,F)+C{M) f \\Kt\\min) + \\^lv^\\minJds. (9.7) 

Jo '- ^ 

By choosing 6 = 5i = 82 = 1/8 and T > so that CT^/'^ < 1/8 in dH]), we find 
that 



t ^ 



I, 



\\vAm + NkIIh^c^) ds < CN3{uo,F) + -\\V^oh.\\h(^r) 



' '|V2/z,,ll2 



ds. (9.8) 



8' 
+ C{M) / ||V^w„|j^i(j^^) +-f'^(s)(||V;j/i„||i2(r) -r II vo»«tiiL2(r)^ 

Combining the estimates fTTf^i . ((8J)) . ([O)) and (H^]) with ([QJ]) . 

II"k|Ihi(o) + ll^o^KllL2(fii) + ll^o^Kllff^fr) + 1 1 ^Kt 1 1^2(0) + ||Vo/iKt||i2(r) (i) 



l|Vi^K||i2(fj) + ||VVot;K||i2(o,) + W^^WWhi^n,) + \Kt\\m 



m 



ds 



<C'N3iuo,F) + C" J [\\v.t\\l^n)+Kis)(^\\VtK\\l.^r) + \\'^'oh^t\\h^T) 
for some constant C" depending on M and e. By the Gronwall inequality and 



ds 



sup 

0<t<T 



l^'«lli/2(o) + WvKtWhift) + l|Vo/iKt|li2(r) + IIVq/ikII 



L2(r) 



+ IkKllHi(n) W + ll«..|lv.^(T) + ll'7,.|li2(o,T;H2(o)) < Cie)N3iuo,F). 
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9.3. Weak limits as k ^ 0. Just as in Section [8.81 the weak limit {v^,h^,q^) 
of {v^,h^,q^) as K -> exists in V{T) x L'^{0,T-H^[T)) x L^{Q,T;H^{n)) with 
estimate 



sup 

0<t<T 



+ hc{t)\\min)] + II«kIIv3(t) + ll9.|li2(o,T;ff=(o)) < C{€)N:i{uo,F). (9.9) 
dnH) imphes that for a.a. t e [0,T], 

IkJOIU-cn < C'(i) 

for some C{t) independent of k, and therefore for a.a. t G [0,T], 

K / Aq^k • '^ofdS -^ 
as K -^ 0. This observation with (|8.20p shows that {ve,hf,q^) satisfies, for a.a. 

te[o,r], 






h^„ o yy'^iyj'^ + (y9^ 



dS* 



(9.10) 

for all iy9 € iJ^'^(ri; F). Since (|9.10p also defines a linear functional on H^{il), by the 
density argument, we have that (|9.10p holds for all (p E H^{il). As {ve,h^,qg) are 
smooth enough, we can integrate by parts and find that {vg,h^,q^) satisfies (|7.2p 
with (|7.2cp replaced by 

[i^D^{v,)i - qeSi]a'^Ni = a[e[(£,- (/le) +M{h)){VoK -1)] ° V^] on (0, T) x T . 

(9.11) 

9.4. _ff^-^-regularity of /i„. By ()9.1ip . we have 
Lemma 9.1. For a.a. t e [0,r], h,{t) e H^-'^iT) with 

||/i.|lH5..(r) < C{M) [||«.J^i(f,) + ||Vi;e||i2(o) + W^lvefm^n,) + llV^/i.Hi^fr) 

1 



and hence 



+ II^IIhmo) 



1^ l|2 



<C{M)e^^^'^+^N3iuo,F). 



(9.12) 
(9.13) 



Proof. We write the boundary condition (|9.1ip as 

CniK) = ^J--\-Woh, 1) • {e-i [[j^Ai(«.)i - g.<J^]a,'A^£] } o r' - -Ml/*)- (9.14) 

By Corollary 17.11 £;- is uniformly elliptic with the elliptic constant vi which is 
independent of M which defines our convex subset Ct{M). Since h G ?i(T), 
M{h) G L^{0,T;H^-^{T))r]L°°{0,T;H\T)), and hence by (15TB . the right-hand 
side of (|9.14p is bounded in iJ^-^(r). The important point is that these bounds are 
independent of e. Thus, elliptic regularity of Cj^ proves the estimate 



||/^e||l,5.5(r) < C{M)[\\D^{v,)\\%^,.^^) + WqeWni-.^D + 1 
so that with ([8^ . ([9?T2l) is proved. 



D 
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9.5. Energy estimates which are independent of e. Having estimate (|9.12p . 
one can follow exactly the same procedure as in Section 19.21 to show that the 
constant C appearing in (|9.9p is independent of e, provided that we have an e- 
independent version of (|9.4p . By Appendix IB] 2. we indeed have such an estimate: 



'^i 






2 II ■ u-t\-/iii.-(r) 
+ CN2{uo,F) + C I K{s)\\V''oK\\l2(r)ds + iS + Ct^/^) I \\v,fHsr^.ds 



^ ds + {S + Ct'/^) f ||«.||l,3(o)^ 
Jo 

Jo 
for some constant C depending on M, S and 6i. Therefore, we can conclude that 

\h(n) + W^etWh^n) + l|Vg/ietl!i2(r) + llV^/i.Hi^fr) (9.15) 

+ Ikell^Hf^)! W + ll"ellv3(T) + he\\h^o,T;HHn)) < C{M)e^^''^+^ m{uo, F) . 



sup 

0<t<T 



Remark 15. Literally speaking, we cannot use Vq(Ci VqW^) as a test function in 
119.1 (A) since it is not a function in H^{fl). However, since /ig £ H^'^(T) for a. a. 
t e [0,T], (MW also holds for all tp G H^{ny (1 H-^-^{T) and Vg(CiVgue) is a 
function of this kind. 

9.6. Weak limits as e —> 0. The same argument leads to that weak limits of 
{ve, he,qe) (dcnotcd by {v, h, q)) as e ^ exists and {v, h, q) satisfies (|7.1I) . 



9.7. Uniqueness. In this section, we show that for a given {v, h) G Yt, the solution 
to (|7.ip is unique in Yt. Suppose {vi,hi) and {v2, h^) are two solutions (in Yt) to 
(|7.3p . Let w — vi — V2 and g = /ii — /12, then w and g satisfy 



{wt, <^) + 2 / Df,w : Dfjifdx + a Lj( {h,aWa - Wz)ds\ 



O T] X 



X {~-h^c.ofj^ip°' + ip^)dS^O 



(9.16) 



for all ip e Vv{T) with w{0) = 0, where L equals L except Li = L2 = 0. Since w is 
in Vv{T), letting w — (p \n ()9.16p leads to 



\\vfH^^^n) + \\'^U\l^n,) + \\<h\\ 



L^T) 



ht\\l2^n) + \\'^lht\\l-i,r) (t) 



t ^ 



l|Vf||i2(o) + !l VVoz;||i2(^^) + ||VV2i;|||,(f,^) + WvtWli^n) 



ds 



<C{M) K{s) ||V^/i||i.(r) + ||V^/jJi.(p) 



ds. 



Therefore, by the Gronwall inequality and the zero initial condition (w(0) = 0), 
we have that w (and hence g) is identical to zero. 

10. Fixed-Point argument 

^From previous sections, we establish a map 0t from Yt into Yt, i.e., given 
{v,h) S Ct{M), there exists a unique QT{v,h) = {v,h) satisfying (|7.ip . Theorem 
14. H is then proved if this mapping 8t has a fixed point. We shall make use of the 
Tychonoff Fixed-Point Theorem which states as follows: 
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Theorem 10.1. For a reflexive Banach space X, and C d X a closed, convex, 
bounded subset, if F : C -^ C is weakly sequentially continuous into X, then F has 
at least one fixed-point. 

In order to apply the TychonofF Fixed-Point Theorem, we need to show that 
Q{v,h) e Ct{M) and this is the case if T is small enough. In the following dis- 
cussion, we will always assume T is smaller than a fixed constant (for example, say 
T < 1) so that the right-hand side of (|9.15p can be written as C{M)Nz{uii, F). 

Remark 16. The space Yt is not reflexive. We will treat Ct{M) as a convex subset 
of Xt and applied the Tychonoff Fixed-Point Theorem on the space Xt ■ 



Before proceeding the fixed-point proof, we note that lemma [6731 implies that for 
a short time, the constant C{M) in (|8.1[) and (|8.4p can be chosen to be independent 
of M. To be more precise, for almost all < i < Ti, 



|9|Il2(s:2) 



<C 



\vt\\l2(n) 



l^vWhm 



Wlhf 



L^iT) 



WWh-^iii) 



(10.1) 



and 



hfuHa) + hfuHu) < C WvtWl.^n) + llVf|l?,i(o) + W^^Amin,) (10.2) 



+ l|V^«|||,i(a) 



2 < r 



-ii^ir 



for some constant C independent of M . 



V7 l|2 

|Vw||l2(0) 



1 



W\\H^{il.) 



Y72 ||2 



\<hr 



L2(V) 



(10.3) 



10.1. Continuity in time of h. By the evolution equation (j7.1dp and the fact 
that V e V3(Ti), ht e L^{Q,Ti;H^-'{T)). Since h e L^{0,Ti;H^-^{T)), we have 
that h £ C*'([0,Ti]; i?'*(r)) by standard interpolation theorem. Although there is 
no uniform rate that h converges to zero in iJ^(r), we have the following. 

Lemma 10.1. Let{v,h) = QTi{v,h). T/ien ||/i(t)||j:/2.5(p\ converges to zero as t ^> 0, 
uniformly for all {v, h) G Cti (M). 

Proof. By the evolution equation (|7.1dp . 

\\hmH2-^r) < I \\h,c.Va~Vz\\m-^o(r)dS<C{M)N3{u„,F)^/V/^. 
Jo 

The lemma follows directly from the inequality. D 

By lemma 110.11 and the interpolation inequality, we also have 

Lemma 10.2. ||Vg/i(<)||//i.5(p) converges to zero as i ^ 0, uniformly for all h G 
Cti (M) with estimate 



\Wlh{t)\\Hi-^r) < C{M)N3{uo,F)t'/^ 



(10.4) 



for allO<t<Ti. 
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10.2. Improved energy estimates. In order to apply the fixed-point theorem, 
we have to use the fact that the forcing F is in V^{T). We also define a new constant 

N{uo,F) := ||uo|ll,2.5(^) + ||F||?,2(^,) + ||^|li=.(o,T,;L^(n)) + ll^(0)ll?,i(a) + 1- 
Note that N3{uo,F) < N{uo,F). 

Remark 17. For the linearized problem ^7.1\ ), we only need F e V^{T) to obtain 
a unique solution {v,h) G Yt ■ 

10.2.1. Estimates for "^"f^v near the boundary. Note that 

1 d 

2di 



CiVgi.||i.(o)+fT/ QBA^P'''^lh^p^lh,,sdS +^||CiA-,(V^«)lli2(,j) 



{F^VliCfvlv))^'^ \vlidtd^y, + Vlia'l~a^)v:Mc!ylv^),kdx 



n n 3 8 



fc = l k = \ 



where /fe's and Jfe's are defined in Appendix [Bll (with ~ replaced by ", and no e 
and ei). 

As in [TU] and [TT], we study the time integral of the right-hand side of the identity 
above in order to prove the validity of the requirement of applying Tychonoff Fixed- 
Point Theorem. By interpolation and (|9.9p . 



Jn 

t 



<C I ||aa||^2(j^)||Vw||ioo(o)||w||ff3(a)ds 
Jo 

< C{M)C{5)N{uo,Ff'^j WvW'-ji^.^^ds + SCiNmiu^^F) 

< C{M)N{uo,F)\c{d)t^/^+6 . 



Similarly, 



Jn 

t 



"'SI 



Vo(afaf)Vowf, + Vo(5fa^^)Vo<,] {(:lylv^),kdxds 
D.f,{Vlv)la'tCiCi,kVlv'dxds < C{M)Niuo,F) \t^'^ + C{5)t + 5 
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For the pressure term, by interpolation and (|8.10p . 

Jn 
< C{M) |k||L=°(o) + Ikllwi.i(fi) + lkll_ffi(o) \\v\\m{n)ds 



<C{M)C{6) / \\q\\HHn)ds + S Mlho,T;H^Q) + h\\L^o,T;m(n)) 
Jo '- 

< C{M)N{uo, F) [c((5)ti/2 ^ g 

By the estimates aheady estabhshed in Appendix iBl with the help of (|6.6p . it is 
also easy to see that 

t 3 8 

{Y,Ik+J2 "^'^)'^* ^ C{M)N{uo, F) [ii/4 + f/2 ^ c{S)t^/^ + 6 

k=l k=l 

Finally, for the forcing term, by the extra regularity we assume on F, 

ft ft ft 

{F,Wl{i:fWlv))ds < / \\F\\H2^n)MH-(n)ds<N{uo,F)+ \\v\\jj,^^)ds 
Jo Jo 

< N{uo,F) + C{M)N{uo,F)t. 

Therefore, 

'- ^ Jo 

< II"o|Ih2(j^) + CN{uo, F) + CiM)N{u„,F) [c(^)(i3/4 ^ ^2/3 ^ ^1/2 + i) + ^ 

By Corollary O 



l|Vg«(i)||i.(n,) + l|V^Mi)lli2(r 



\'^ov\\m{n,)ds 



<CN{uo,F)+C{M)N{uo,F) C{S)0{t)+S as t^O (10.5) 

where C depends on v, a, vi and the geometry of F. 

By similar computations, we can also conclude (the ()7.27|) . ()8.9|) and ()9.5|) vari- 
ants) that 



Ht)\\l-(n) + \\^lKml^(T) 



+ / hfm{n)ds 



<CN{uo,F) + C{M)N{uQ,F)0{t) as i -^ ; 



l|Vo«(i)|li2(o,) + l|V;]/^(t)li2(r) +/ ||Vo«||^i(^^)ds 

-' Jo 

< CN{uo,F) + C{M)N{un,F)0{t) as i^O; 



l!Vv(i)||i2(o) + / \\vt\\i2^n)ds 
Jo 

< CN{uo, F) + C{M)N{uo,F)0{t) as i ^ 
where C depends on v, a, vi and the geometry of F. 



(10.6) 



(10.7) 



flO. 
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10.2.2. L^H^-estimate for vt- For the time-differentiated problem, we are not able 
to use estimates such as those in sections 18.61 and 110.2.11 since no e-regularization 
is present; nevertheless, we can obtain estimates at the e-regularization level and 
then pass e to the limit once the estimate is found to be e- independent. We have 
that 



^ %t\\Un) + '^\\D,v,\\l.^,,) + ^^^ I QA"P-^'h,^^ph,,,sdS 



2dt" ""'^^"> 2" " '"^^^"^ 2dtjr 

{Ft,vt)-iy / id'la'j)tv^j, + {a'ld'ij)tvli v^j^dx + / qtd{^v'^^dx 






^A^t(^){A''^'^%h,^p 



,7(5 



httdS 



Q^^A'^P^^t.^phu^sdS - / Q^^sA"''^^ht,o,phudS 



e 



-a/37; 



L^P^h^^p^ httdS- / &{L2)thudS + Ki+K3 + IU + K5 + Ke 



where K'^s are defined in Appendix [Cl (withoutei ) . 

As in the previous section, the time integral of the right-hand side of the identity 
above is studied. It is easy to see that 



{Ft,vt) - i/((afa^)tw;'^ -|- (a*'a^)tw>) Wj^fe + Ki + K5 + Kq 



ds 



< CiM)N{uo, F) t^'^ + ^1/2 ^ c{5){t^/'^ +t) + 5 
and by Appendix [Cl particularly Remark 1211 



JT 



e 



2^ " '■''^ '■"' v/d^^ 



v/ditM(^"''^')t^,«/3 

dSds 



,7(5 



< C{M)N{un,F)t^'^. 

Special treatment needed to be done for the rest terms, and we break this procedure 
into several steps. 

Step 1. Let Bi — I I {qdl.)tv^jdxds. By the "divergence free" condition (j7.2bp . 
Jo Jn 



Bi — I I aj^f.qVf gdxds — / / af.^qtv '^dxds. 



Jn 



Jn 



By interpolation and 



akt^v^pdxds 



Jo. 



< C{M)C{5) / |k|li2(^)ds -I- 5 |k|li2(o,T;Hl(n)) + ll"tlli2(0,T;Hla2)) 



< C{M)N{uo,F) C{5)t + 5 . 
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For the second integral, we have the following identity: 



aitQtv'^idxds = / {aitqv';g){t)dx - / 4t(0)q(0)4,fda; 
Jn Jn Jn 

t 



i^kt^.LStldxds. 



Jn 



By the identity a^^ = -~a],v^aj 



(aj,(W f)tqdxds 



JVl 



< 



Jn 



aittv% + di,vlt 



dxds 



<C{M) / (l + ||St||Hi(o))||Vw|U4(o)||<z||L4(o)rfs. 

JQ 



Therefore, 



{ait'v^i)tqdxds 
/o Jn 

<C{M)C{6)N{uo,F) 

< CiM)N{uo,Fy\ciS)it + t'^)+6 



* hrHHn)h\\%';^ds + S f\l + \MJ,.^,,^)ds 

"'O 



3 1 

where a ^ - if n = 3 and a — - ii n — 2. 

4 2 



The second integral equals / Vuq : (Vuo)"^'?(0)dx which is bounded by CA^(mo, -F). 

Jn 

It remains to estimate the first integral. By adding and subtracting / ai^f.{Q)qv '^dx, 

Jn 
we find, by at{0) G H^i^l), that 



/ i4t'iy':i)it)dx 
Jn 



< 



dx 



{di-aimiiv'im 

< C\\dt{t) -at(0)||L4(o)||q||L2(o)l|Vw||L4(Q) 
+ C((5i)||Vi;||i.(j,)+<5i||g||i.(j,). 



aiMl^j 



dx 



Noting that 



\Vv\ 



L2(0) 



= IJVuo 



yvtdsW 



L2(n) 



< 



* -,2 

||Vuo||l2(!:2) + / ||V-Uf||L2(j^)ds 



< 2 



l"o||l,i(,,)+C(Af)iVK,i^)i 



dnH), dnSt) and pini) Imply 

{aitqv\{t)dx < C{M)N{uo, F)t^/^ + C((5i)7V(wo, ^) 



5i 



kt|lL2(f^) + ||Vo/i||i2(r) 



Summing all the estimates above, we find that 

\Bi\ < C{Si)N{un, F) + C{M)N{un, Ff \c{5){t + t'-^)+5 



Si 



\vt\\i^n) + \mh\\hir 
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Remark 18. It may be tempting to use an interpolation inequality to show that 
q £ C{[0,T]; X) for some Banach space X by analyzing qt via Laplace's equation. 
The problem, however, is that the boundary condition for qt has low regularity 
L2(0,T;i/-i-5(r)) (by the fact that ht G L'^[Q,T; H'^\T))), and thus standard 
elliptic estimates do not provide the desired conclusion that qt G L^(0, T; iJ^(f2)') 
(and hence by interpolation, q G C([0, T]; -ff'''^(f2)). However, suppose that qt G 

L^{Q,T; H^{iiy); then we can estimate / / af.tqtv\dxds by the following method: 

Jo Jn 



fl 



cLkt^tv idxds 



< 



\"'W,ia^jv'^i\\HHn)ht\\HHnyds 



< C{M)N{uQ,F) 



t + f 



5/8 



Step 2. Let ^2 = / / 9 

'0 Jr 



[L'^ h,ai3-i\thtt + {L2)thtt dSds. It is easy to see that 



Jr 



e{L2)thttdSds < C{M) / Ml^^t) + ^tlk^cr) 



t ^ 



ds 



<C{M)N{uQ,FY'^{t + t'^/^). 



For parts involving Li, we have 



e 



Jr 



[L'^'^^h^ai3'i\thttdSds = / / 6 



- a/37 



Jr 



Li h^af3yhttdSds (= ^2) 



QL"'^'' ht,ai3jhttdSds. (= B2 



By interpolation, 



\Bl\<C(M) / \\Q\\L^{r)\\h\\w^A^r)\\htt\\LHr)dSds 
Jo 

ds 



<C{M) \\v\\m(^n) + \\vt\\m(^n) 
Jo '- 

< C(M)iV(uo,F)i/2ti/2 
while by ()6.6p and Corollarv l6.H 

\Bl\ < / ||0|ki-5(r)ll^t||H2.5(r)||Li''''|lHi.5(r)||/itt|lffo.5(r)ds 
Jo 

< C(M)||L"''''||^i.5(r) / ||/i||H2.5(r) \\v\\m(n) + \\vt\\m{n) 

Jo '- 

< C{M)N{uo,F)t^/^. 



ds 



Therefore, 



Step 3. Let B3 = / K^ds = 

Jo 
and L2 part of i?3 is bounded by 



IB2I < C{M)N{uo, F){t + ^3/4 + ii/4). 



C{M) 



fmiH 

Jo 



e[L-^ih)]t[iv o ,7--) . {Voht)]dSds. The L^ 
i-^(r)l|w||i/i-5(r)l|/i|k3-5(r)ll^t||i/2(r)||^tlk2(n)rfs 
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and hence 



e 







L'^^'^h 



aPj 



L, 



[{v o f]-^) ■ (Voht)]dSds < C{M)N{uo, F)t^/^. 



By the integration by parts formula, the highest order part of B-^ can be expressed 
as 



* <" e{vofj-'') 



Jo Jr 



\/det(5o) 



V^drt(^(A"'^''^)t/i,„^ W^htdSds {=Bl) 

J ,7(5 



+ / Q{vof]-^)A"''^^ht.,c.pVoht,^sdSds {= bI) 



[Q{v o r')],^A'^P^^ht,o.pVoht,5dSds (= Bl) 



JT 
t 



JT 



'0 Jr 

It is easy to see that 



Mv o Tl].nsA'^'^''^Kc.ti^^htdSds. (= Bt) 



B'^\ < C{M) / \\evof^-\\Hi.^r)\\ht\\HHr)\\h\\HHr)\\ht\\HHr)dS 
Jo 
< C{M)N{uo,F)t 



and 



\Bl\ < C{M) / \\Qv o f]-^\\wiA(T)\\A\\L^(T)\\ht\\m(r)\\ht\\w^A(r)dS 
Jo 

For _B|, by the integration by parts formula, 



Bi 



1 

2 7o 



e(wo^-^)i"'3T"5Vo 



ht.apht.-yS 



2 Jo Jr v/det(go) 



dSds 



ht^apht,fsdSds 



and hence 



\Bl\ < 



l|Voe||L4(r)||wA||i^(r) + ||6||L~(r)l|w-4||;yi,4(r) 
X \\ht\\w2A(r)\\ht\\H^{r)ds 

<C{M)N{uo,Fy/^ f \\v\\H^n)ds 
Jo 

< C{M)N{uo,F)t^^^. 



For B^, noting that 



9,^5 = det(Vo?7^),75v/det(G,-) o 77^ + det(Vo^^),-y Vdet(G^) o ^^^^ 
+ det{Vovl.sV^^M(h>¥,^ + det(Vo^")Vdet(G,-)oJ7-^^^ 
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and ||Vodet(Vo^^)||H0.5(r) < C{M)t^/'^, wc find that 

\Bt\ <C{M) / ||Vodet(Vo7?")||H0.5(r)||V^/it||H0.5(r)||Vo/it||//i-5(r)ds 

+ C{M) / II det(Vo?7")l|L~(r)||Vo?7"|li~(r)l|Vo/itlU^(r)||Vo/it|U^(r)ds 

"'0 

< C{M)N{uo,F)t^/^+C{M)N{u^,Ff'^ f \\v\\]is^^)ds 
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Combining all the estimates, we find that 



ft .t 

Step 4. Let S4 = / K^ds = / / 9 



grating by parts, 



JT 



Lj,{h) [{Woh, -l)t ■ {v o fj~'')]dSds. Inte- 



Bi=- I I Lnih) [Qti^oh, ~l)t ■ {v o fj--) + &{\/oh,-l)t ■ {v o f]--)t 

+ e(Vo/i, -i)« • {v o 7;"")] dsds + I eL^(/i)[(Vo/i, -i)t • {v o f]-^)]ds. 

For the first integral, (|6.8p implies 

j QL~^{h)[{Voh,~l)f{vorn]dS 

< ||e||L-(r)||i,;(/^)||L2(r)||Vo/it||L4(r)||uo7?"^||i4(p) 

<C{M)N{uo,F)\\ht\\H^..^r) 

< C{M)N{uo,F)t^/^. 

It is also easy to see that 



Jr 



L-hih) 



Qt{Vji,-l)t ■ (v o f^--) + Q{Voh, -1), • {v o fj-^\ 



dSds 



< C{M) ||w||L~(r) + ||«t|U*(r) ||-^,;(^)||L2(r)||Vo/it||Li(r)c's 



< C{M)N{uo,Fy/^ / [||z;||«3(^) + \\vt\\H^in) 



< C{M)N{uo,F)t^/^. 
For the remaining terms, H'^-^(r)-H~'^-^{r) duality pairing leads to 

/ j QL~^{h){Voh,-l)tfvdSds 
Jo Jr 

< 



\Q\\mHr)\\L~^{h)\\HO-^{r)\Mm-Hr)\\htt\\Ho^r)ds. 



By interpolation. 



1/2 



1/2 



\Lj^ih)\\HO-^r) < C{M) ||/^||;;.5(r)||/i|i;;,^5(r) + 
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and hence 



"T 



eL~Jh)(Woh, ~l)tt ■ {v o 7^-^)dSds 



ds 



< C{M)N{uo,F)J \\htt\\HO-^ir)[\\^lh\\Y^^^^ + 

<C{M)C{S)N{uo,F) f \\\\/lh\\L2^r) + '^]ds + 6C{M)N{uo,F) 
Jo '- 

< C{M)N{uo, F) \c{S){t^/^ +t) + S 
All the inequalities above give us 

\Bi\ < C{M)N{uo, F) [c{S){t^/^ +t)+ ii/8 + s 
Summing all the estimates above, we find that 

Jr ^ Jo 

< ll«*(0)|li.(o) +^l \GfhtMO)\^dS +iC + Ci6^))Niuo, F) 

+ C{M)N{ua, F) \c{d){t + t"^'^ + t^'^ + t^'^ + t^/"^ +t^) + S 

+ Si[\\vt\\l2t^a) + \\^th\\hir) 
and by Corollarv l7.ll 

Ivtimhin) + llVg/»t(i)lli2(r) 



ds 



\vt\\H^n)ds 



< [C + C{Si))N{uo, F) + CiM)Niuo, F) CiS)Oit) + 6 



(10.9) 



<5i 



l^tllL2(r2) + W^oHl^if) 



where C depends on ly, a, vi and the geometry of T. Since this estimate is inde- 
pendent of e, we pass e to zero and conclude that the solution {v, h) to (|7.ip also 
satisfies ((TU^. 

10.3. Mapping from Ct{M) into Ct{M). In this section, we are going to choose 

M so that e{v,h)e Ct{M) if {v,h) c Ct{M). 

Summing pUS]) . (fTUTB)) . pHT)) . (fTUTS]) and (HHH), by dES]) we find that 






t ^ 



hfm(n) + \\^ov\\]ii(^m) + l|Vow||l^i(Oi) + htW'H^in) 



ds 



<{C + C(<5i))iV(wo, F) + CiM)N{uo, F) C{S)0{t) + 6 



■Si 



IvtWhi^a) + 11^0^111,2(1 



NAVIER-STOKES INTERACTING WITH A FLUID SHELL 



43 



where C depends on u^ a, vi and the geometry of F. Choose 5i 



2' 



<tWmn) + l|Vo^'W||i2(o,) + l|Vg«(i)||i2(o,) + \\vt{tWmn) 
+ l|Vg/^(t)||i.(r) + ||Vg/i(i)||i.(p) + ||V^/j(i)||i.(r) + ||V^/i*(t)||i.(p) 



+ 



\Hm(n) + I|Vo«|Ihi(,,,) + ||V2z;||2^i(j,^) + ||z;t||2,,(^^ 



ds 



< CiN{ua, F) + C{M)N{uo, Ff C{5)0{t) + 5 



where Ci depends on v, a fi and the geometry of T. Shnilar to Section 18.71 for 
almost all < i < T, 



Hmlnn) + \\vt{t)\\l2^n) + W^lHrnlnr) + l|V^/^t(i)lli2(r) 



lkllff3(jj) + ||ut|lHi(o) + lkllH2(s:2) 



ds 



(10.10) 



< C2N{uo, F) + C{M)N{uo, Ff C{5)0{t) + S 



for some constant C2 depending on Ci 
By dUD, ^^ and (FTTdl) . 

\\ht\\']j2-5iT\ds < 



5(r)^ 



l + ||/i||^3.5(r) ||w|ll/2.5(r)rfs 



<C{M)N{uo,F)t^/'^ 



(10.11) 



and 



\\htt\\%o.5^r)ds < C{M) / ||/it||l/i-5(r)llw||H2(s:2) + l|/i|ll/2.5(r)||wt||^i(o) 
"'0 '- 



<CiM)N{uo,F) 
Also, by (fTU:^ and (fTinUl) . 

|2 



tl/4+^1/2 



ds 
(10.12) 



||/i||^5.5(r)C^s < C 
"'0 



vt\\m(n) + I|Vw||l2(o) + l|Vov|lHi(ni) 



iv^/^r 



l^ll 



1 



ds 



l//i(0) 

< C3^(uo, i^) + C{M)N{uq, Ff C{5)0{t) + 5 



L^r) 



(10.13) 



for some constant C3 depending on C2. 

Combining (jlO.lOp . (jlO.Up . (|10.12|) and (|10.13p . we have the following inequality: 



HtWuHn) + html2(n) + ||Mi)llH4(r) + \\ht{t)fH2(r) 



II«IIh3(o) + htfrn^i^) + ||/i|lH5.5(r) + ||/it|lH2.=(r) + \\htt\fH«-'-(v) 



< {C2 + C3)N{u„,F) + C{M)N{u„, Ff C{6)0{t) + S 



ds 



Let M = 2(C2 + C3)N{uq, F) + 1 (and hence corresponding Tq and T in Lemma 
and Corollary [73] are fixed). Choose S > small enough (but fixed one such 5) 
so that 

C{M)N{uo,Ff5< - 
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and then choose T > small enough so that 

C{M)N{uo,FfC{6)T<]. 



Then for almost aU < t < T, 



l«WllH^(i2) + ll«tWI" 



L^{n) + ll'^(OII_H-4(r) -T ii"tv'-^iiH2(r) 



\ht{t)r 



lkll_ff3(o) + ll^tll_ffi(s^) + \\^t\\H^-^r) + \\htt\\HO'-{r) 



ds 



<C2Niuo,F) + - 



and therefore 
sup 

0<t<T 



\v{t)f 



\\vt{t)f 



m(n) -t- w^ty-jWL-^^n) + \\Ki)\\m{T) 



\ht{t)f 



m{T) 



+ lkllv3(T) + ll/^ll^^m < 2C2A^(uo,F) + 1, (10.14) 

or in other words, 

Remark 19. jlO-lj^ implies that for {v,h) e Ct{M) (with M and T chosen as 
above), the corresponding solution to the linear problem |7.i[ ) {v,h) ~ 0t(w,/i) is 
also in Ct{M). 

10.4. Weak continuity of the mapping 8^. 

Lemma 10.3. The mapping Qt is weakly sequentially continuous from Ct{M) into 
Ct{M) (endowed with the norm of Xt)- 

Proof. Let {vp, /ip)pgN be a given sequence of elements of Ct{M) weakly convergent 
(in Yt) toward a given element (w,/i) G Ct{M) {Ct{M) is sequentially weakly 
closed as a closed convex set) and let (^^(p), ^CT(p))peN be any subsequence of this 
sequence. 

Since V^{T) is compactly embedded into L^{Q, T; _ff^(J7)), we deduce the follow- 
ing strong convergence results in L^{0, T; L'^{il)) as p ^ oo: 

{aj)p{a^)p -^ a\a\ and {a'i)p{ai)p -> a\ai, (10.15a) 

[(ai)p(4)p]j ^ (a£4)j and [{al)p{ai)p]^j -^ (a^aDj > (10.15b) 

{a'Dp ^ al (10.15c) 

Now, let {wp,gp) = Qrivp^hp) and let qp be the associated pressure, so that 
('Zp)peN is in a bounded set of V^(T). Since Xt is a reflexive Hilbert space, let 
(u'o-(p),5o-(p)5 9ct(p))p6N be a subsequence weakly converging in Xt x V^(T) toward 
an element {w,g,q) £ Xt x V^(T). Since Ct{M) is weakly closed in Xt, we also 
have {w,g) G CTiM). 

For each G L'^{0,T;H^{n)), we deduce from ((7?^ (and RemarklH]) that 

{wt,(l))L^(n) + ij DriW : D,j(j)dx + o" / Lh{g){g^a4'a - 4>z)dS 



qal(j)'^jdx 



dt 



{F,^)dt 



which with the fact that, from (|10.15|) . for all t G [0,T], w G Vi,, provides that 
{w,g) is a solution of (|2.16p in Ct{M), i.e., {w,g) = QT{v,h). 
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Therefore, we deduce that the whole sequence (0T(wn, hn))neti weakly converges 
in Ct{M) toward 8t(w,/i), which concludes the lemma. D 

10.5. Uniqueness. For the uniqueness result, we assume that uq, F and T are 
smooth enough (e.g. uq G H^-^{n), F G V''(T), T is a H^'^ surface) so that 
Uq and the associated Ui, q^ satisfy compatibility condition (|4.4p . Therefore, 
the solution (u,/i,g) are such that v G V^{T), q G L'^{Q,T;H^{n)) and h G 
X°°(0,r;ij7(r))nL2(0,r;7j8-^(r)), ht G L°°{Q,T]H^{T))f^L'^{Q,T■,H'^■^{T)), hu G 
L°°(0,r;ij2(r)) nL2(o,T;i73.5(-p)) This implies a G L°°{0,T; H^{n)) and hence 
by studying the elliptic equation 



{aia'iqt,k).e = K(apa^^'j),M + 4^/ + ^i^/ " [(ai<)t9,fc],£ 



in n. 



g, = J-' CTLH{h)N, - z.i?,(z;)faf TV, - K^V,)*? 



afTVf on T, 



we find that qt G L^iO,T; H^{Q)) and this implies «« G L'^{0,T; H^{Q)). By the 
interpolation theorem, we also conclude that Vt G C'^([0,T]; iJ^'^(ri)). 
Suppose (w, ft., q) and (u, /i, g) are two set of solutions of (|l.ip . Then 

(w - v)t - iy[4D^{v ~ w)^],fc = - af (g - g),fe + SF 



al{v — vYj — 6a 



iy[D^iv-v)]i^iq-q)S: 



a\N, 



uQ 



V 



L„(/i-/i)(-Vo/i,l) 
+ 5Li + 5L2 + SL3 

{h - h)t O if = [h,a O 77'"](Wq - Va) ~ {Vz " V^) 

+ 5hi + 5h2 + (5/i3 
(w-w)(0) = 
(/i-/i)(0) = 



where 



dF = fori- fofi + v[{a'la] ~ d!}hl)v'j\,k + v[{a'la{ - a^a^){)^,],fc 
5Li = aQ 



5a ^ {al-hly. 



LhCh)(^oh - Woh, 0)J o if - v[a\a\ - a^hj)v%Nj 
- ^(4fl£ - 4ai,W,kN, + (ai ~ ai)qNj 
SL2 = Q[Lj^{h) o ?7"](Voft o ri^ - \/oh o 77^ 0) 

BLhCh) o Tj^ - BLhCh) o fj^] (Voft o fi\ -1) 

SL3= \[Lhih)-L'^{h)]{VoK-l) 
Shi — {h,a o rf ~ h^a ° ff)Va 



o ff 



Sh2 = 



{h,a - ft, a) O??^ 



(5ft3 = - {ht o 7]'^ - ht o fj'^) 



3 3 

We will also use SL and 6h to denote J2 ^k and J2 ^^k respectively. 

k=l k=l 



(10.16a) 
(10.16b) 
(10.16c) 

(10.16d) 

(10.16e) 
(10.16f) 

(10.17a) 

(10.17b) 
(10.17c) 

(10.17d) 

(10.17e) 
(10.17f) 
(10.17g) 
(10.17h) 
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Similar to (11.3) in [TT], we also have the following estimates. 
Lemma 10.4. For / e H^{fl) and g £ H^'-'iT), 

\\f on- f o fj\\L2m) < CVi\\f\\H2(n) / \\v - v\ 



\v - v\r„n^^ds 



llgorf -gor7||i2(r) < CVt||g|Ui.5(r) / \\v - v\\hi (^n)ds 
for some constant C. 



1/2 

7 

1/2 



(10.18) 
(10.19) 



Remark 20. Assuming the regularity of h, hf and hu given in the beginning of this 
section, we have 



ll^-^2||//2(r) + ||(5/ii +^/i3||//2.5(r) < C\/t / ||w- w||^3(o)rfs 



1/2 



and 



\\{5L2)t\\mT) + \\{5hi+5h:i)t\\H^(j) 



< c 



\v-v\\HHn) + 



V~t 



"''li/2(Q)l 



1/2- 



(10.20) 
(10.21) 



\v~v\\H^(n) + \\v-v\\H3{n) 



+ \/i||/itt||ff3-5(r)( / h - v\\H^n)dsj 
By using (|10.18p to estimate ||5i^||i2(j2), we find that 



1/2 



\\yiv~vm\\hin)+ I Uv~i)t\ 



2 ^« 



< C{5) / \\v - v\\l,^n. + \\h- hWJ,,.^. ds + {C{5)t^ +5) \\v-i\ 



(10.22) 



Lff2(ri) 



ds 



Uv-v)t\\HHn) + h-i\\m 



(O) 



ds. 



(10.23) 



For the LjH^ estimate for -y — -D and the LfH^ estimate for (v — v)t, we have 



ld_ 
< C 



llCiV^(i; - «)||i2(o) + 2aEniVt,{h - h)) + -\\C,D,Vi{v - i)\\i. 



m 



\\SF\\m(n) + ll(« - ^)t\\i2(n) + l|V(« - v)\\i,(,,^ + \\\/Vo{v - «)||i2(o;) 
+ II V^(/i - h)\\l.,J +S\\v- v\\l,,^. +D,+D2 + D, 



and 



1 d 

2di 

< C 



II (« - ^)t\\Un) + 2'^Eh{{h - h)t) + -||V(« - i)t\\h(n) 



\^t{h - /i)lli2(r) + W^tih - /i)t|li2(r)) + ll-JFiH^.^^,),] + 5\\v - «||^3(a) 
E'l + £'2 + £'3 . 
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where 



Di 



Ci V;^(g - q)Vidadx, D2:= e [Lh{h - h)] o 7?^ {Vl5h)dS, 
n Jr L J 

dL-W*{v-i)dS, 



and 



Ei-.^ {q-q)t{Sa)tdx, £^2 := / e[Lh{h ^ h)] o 7^^ {5h)tdS, 

Es-.^ I {5L)f{v-i)tdS. 
By using (|10.20p to estimate A and (|10.2ip . (|10.22p to estimate Ei, we obtain 



\m{v-v){t)\\h(n,) + \mih-h){t)\\i.^^~^ + / \\yyi{v~v)\\i.^^^)ds 



t ^ 



< C{S) / \\{v- f)),|li.(o) + ||Vo(« - S)||i.(f,) + ||V^(/i - /i)lli.(p) 
+ {C{S)t'+S) f \\v-v\\l.^n)ds + S f \\q-q\\l.^n)ds 



ds 







(10.24) 







and 



iv-v)tmUn) + mih~h)t\\i.^r) + i \\^iv-v)t\ 

„~.l|2 







L2(0)"'* 



< C(J) / II" - nimn) + l|V;5(/^ - ^)lli^(r) + (1 + l|/i«llff-(r)) 







x|iVg(/i-/i)Ji.(p) 



ds 



(10.25) 



+ {C{5){t + t')+S) / |k-S||?,3(o)ds + <5|lg-g|li.(,,) 
Jo 



s 



\\{v-v)t\\m{n) + h-Q\ 



;;l|2 



ds. 



Summing (|10.23p . a0.24p and (|10.25p . we find that 

r(t) + / Z(s)ds < C{S) [ k{s)Y{s)ds + {C{6){t^ +t)+S) [ Z{s)ds (10.26) 
Jo Jo Jo 

where 

Y{t) = [||^; - i{tWH^^n) + ||Vg(« - i)(t)\\hin,) + ll(« " ^)t(i)lli.(o) 

+ ||/i-/i||^4(r) + ||(/i-^)t|ll/2(r) , 
Z{t) = \\{v- iUml^^n) + II VV2(« - i){t)\\l^(n,y 
By letting 5 = 1/4 and choosing r„ < T so that C{d)(T^ + T„) < 1/4, 



Y{t) + f Z{s)ds <C [ k{s)Y{s)ds 
Jo Jo 



(10.27) 



for aU < t < r,(. Since 1^(0) = 0, the uniqueness of the solution follows from that 
Y{t) = for ah < t < r„. 
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Appendix A. Elliptic regularity 
We establish a K-independent elliptic estimate for solutions of 

of,- + K^lv^ = f (A.l) 



e 



\/det(5o) 



/ ,7(5 



where h^ and v^ satisfy dH]) with K & i?*(r) v^ e H'^iT), and / G iJi'5(r). 
Letting w = w^ o ^~^, (|A.1[) is equivalent to 



e 



v/det(5o) 
which implies 

e 



yd^t(^A"^^^/l,,„^ i-Voh, 1) + «A> = / o jy- 

J ,7^ 



V^detCgo) 



v/det(.go)A"'^''^/i«,„^ + kJ^^AIw ■ {-Wo'h, 1) 



,7^ 



(A.2) 



(A.3) 



Recall that w ■ (— Vq/i, 1) = /ik*- 

Let Z?;i denote the difference quotients (w.r.t. the surface coordinate system). 
Taking the inner-product of (|A.3p with D^hDh^Qh,^, by CoroUarv lT.ll we find that 



ds. 



I'll ll^hVo/i«||i2(r)ds < C(e) / ||/iK||H2(r) + ||/||/i-i(r) + «l|w||H4(r) 



Since the right-hand side is independent of difference parameter h, it follows that 
h^ e H^(T) (as it is already a i/^-function) with the estimate 



llV;5/»Jli2(r)ds < C(e) / ||/iK|lH2(r) + ||/|lHi(r) + «ll^llH^(r) ^s- (A.4) 
■/o '- -' 

Next, we obtain a K-independent estimate of K||w|j^4cpi. ^^ taking the inner- 
product of (|A.2p with VqW and VqW, we find that 

73u U\\\2 , „. / ll„„l|2 



l|VS/iK(i)|li2(r)+K / |k)|||,3(r)ds 



< C(6) / ||Vi^/j«||i.(r) + ||/||i.(p) + ||u;||^..5(,,) 



ds. 



and 



|Vo/iKWIlL2(r) +« / lhllH4(r)ds 



(A.5) 



(A.6) 



Jo '- ^ Jo 

where we use (jA.51) to estimate k / ||i(;||/f3/p\(is. (jA.6p provides a K-independent 

Jo 
estimate for KHuill^i.ps; hence by choosing 6i > small enough, (jA.4p implies that 

for aUt e [0,r], 



^lh,\\H-,(r)ds<C' W^th. 



t ^ 



for some constant C" depending on e. 



\LHT) 



fl|2 



1.5(p) 



\H^(n) 



ds (A. 7) 
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Appendix B. Inequalities in the estimates for VqW near the boundary 
B.l. K- independent estimates. Since ^i = 1 on F and 

i-Voh o fi\ 1) • Vlv^ - Vt{{-Vji o f]\ 1) . V,) - VU-^oh o f]\ 1) • v^ 

- 4Sll{-Voh o f]\ 1) . Vot^K - 6V2(-Vo/i o vi\ 1) • V^w, 
-4Vo(-Vo/io77^,l). Vj^i;^, 

wc find that 



e 



L-h{K) o n^ ii-voh o f}\ 1) • vl{c!ylv,))ds 



e 



i/i(^«) ° V^\ [v;5(-Vo/i o 77^ 1) • f;« + 4Vi](-Vo/i o ,7^ 1) • VqWk 



Vr« 



^^(/ik) ° ^"J (Vo(-Vo/i o jy^ 1) . Vlv,)dS (= /a) 

e 



--^t 



Vd^t{^ (Lf-'h^p^ + L2) o fj- VliK, o ,7-)d5 (= 73) 



v/dct(5o) if-'/i^^/s^ + L2 ° ^" V2(/i«, o ,7-)d5 (= 74) 



r Vdct(go) 

2Voe , 

r Vdct(go) 
+ / (Vge) [(Lf-'/i^^^^ + L2 ) o ^^ I Vj^(/i,, o ^-)d5 (= 75) 

e 

r v/dct(.go) 



(yd^tO^A^^-'^/j^,,^),^, o f]^ VliKt o ^-)d5. 



The last term of the identity above, by a change of coordinates, can be written as 
® [( Vd^tMA"'37^/i«,„^),^5 o f]^] VtiKt o f]-)dS 
Vl{^dct{go)A'-P^'K,^p)^^sVlKtdS + Ri 



r v/det(.go) 



r Vdet(go) 
VoG) 



r Vd<3t(.9o 



=Vo 



{./A^t{^)A"P-<'K^^p)^^s o 77^1 v2(/i,, o n-)dS (= Ji) 



( Vd^t(^A"'="'^/i«,„^),^5 o 77" V§(/i«, o ,7-)d5 (= J2) 



/r Vdct(go 
1 rf 

where 73 = 6* (g) 6* (g) 6* (8) 6* with 6 = Vq^^, and 



. y i?A"'5''^v2/7,,„^vg/i«,^,d5 + i?'i 



i?i(t) = / 6* ® 6* ® (V06*) ® (Vo6*)Vo(Vdet(5o)A"'3''^/i«,„^),^5Vo/iKtd5 (= J3) 



b'®b'®h'® (Vo6*)Vo(Vdct(go)A"'3''^/i«,„/3),75V2/i«,d5 (= J4) 
6* (g 6* (g &* (g (Vo6*)V2(Vdet(go)^"''''^/iK,a,3),75Vo/iKtci5 (= J5) 
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and 



-.Wo{Vdet{go)A'^P^')Woh^,c0'^lKt.jsdS (= J7) 



B 



r y/dct{go) 
B 



r \/det(5o) 



S- 



r \/dct(go) 



V^(Vdet(5o)^"''^'>«,a/3V2/i,,^^,d5 (= Jg) 

Vg(yd;tMA"'375/l,,„^)V2/l,,^5d5 (= Jg) 



, ^^^l{./A^t{^)A''P'''K,^p)^lKtdS. (=Jio) 
/r Vdct(5o) 

It follows that 

|/i| < C{e){l + WVtKWL^i^DmWWmin',) , 
\h\ + \h\ + \h\ < C{M){\ + \\h\\HHv)mlv.\\min,) , 
and hence that 

|/i| + I/3I + I/4I + I/5I < C{e) [||V^/i,||2 + ||/^||2^ + ll + S\\v^\\l, 



3(n)- 



It follows that 



IJ2I + IJ3I + IJ5I + iJiol < C{e)\\VtK\\L2^r)\\^lKt\\LHT) 

We need only obtain K-independent estimates for the terms I2, Ji, Ji, J7, Js a-nd 
Jg. By the H^°-^{T)-H°-^{T) duality pairing, 



I/2I < C{M) WWlh^Wn-. 



(r) 



1 



l'i^K|lH2.5(r). 



Therefore, by interpolation and Young's inequality. 



I/2I < C 



''k|Ih4(p) 



1 



'5i||Vo/i«||^3(r) +'5||i'«||^3(a) 



(B.l) 



for some C depending on M, S and Si. 
For Ji, J4 and Jg, we find that 



I Jil + IJ4I + IJ9I < C(e)||/i«,||H4.5(r)||w«|ljy2.5(r) 



<C' 



W^lh^WHHr) + 1 + SiW^lh^Wmir) + S\K\\m 



(O) 



for some constant C depending on M, e, S and 61. 

For J7 and Jg, by the iJ-i-^(r)-iJi-5(r) duality pairing, 

iJr] + \Js\ < C{M)\\B\\Hi^r)\\h\\H^^r)\\hK\\H^^r)\\vK\\H^^r)- 
Similarly to the estimate in (jB.ip . we find that 



|J7| + |J8|<C(M) ||/^.|jl,4(r) + l +Si\\Vlh,\\l^r) + S\K\\m 



(n)- 
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Summing all the estimates and then integrating in time from to t, by Corollary 
17.11 and the fact that B is close to 1 in the uniform norm for T small, 



yllV^/./-"^ 



L2(r) 
t 



-^IJ." 



[Lj,{h^)i~W„h, 1)] o ,7-J . Wt{Ct^iv^)dSds 



Jo Jo '- 

+ ^/ hA^rnds + h I ||Vo/iK||?/3(nds 
Jo Jo 

for some constant C" depending on M, e, S and Si, where 

K{s) :- 1 + ||i)||?,3(o) + ml^D + \\hffH2.^ry 



B.2. e- independent estimates. We next obtain e-independent estimates for the 
first two terms of /i, as well as those for /2, Ji, J2, J3, J4, J5, Jg and Jio with /i^ 
replaced by /ig. Let 



Ii= - I e [Lj^ih,] o fj^^ [VU-'^oh o fj^, 1) • «, 



dS*, 



II 



U^ 



Lj^{h,)ofj- Vli~Voho fj\ l)-Vov, 



dS 



By the iJ-i-^(r)-iJi-5(r) duality pairing, 

|/i| + |/2| < C(M)||L^(/i,)|Ui..(r)|ke|k-(r)|l(Vo/i)o^nk-(r). 
Therefore, by dH]) and ((9l^ . 



|/il + |/il<c(M)ii/4 



'■e|lff5.5(r) 



K'e|l-ff3(^) 



(B.2) 



for some constant C depending on M and S. 

For Ji, we use an L'^-L^-L^ type of Holder's inequality and conclude that 

|Ji|<C(M)ti/2||/i,||H5.5(r)||z;.||ff2.5(r) 

while for the other J terms, we use the i7"-^(r)-i/^"'^(r) duality pairing to obtain 

IJ2I + IJ3I + \Ji\ + \M + \-M + \M < C{M)t'/^\\K\\H^o.r.^r)\\Ve\\H2-Hr) , 

and hence all the J terms are bounded by the same right-hand side of the inequality 
in dEl. Therefore, 



vllVo/^^WI!i^(r) 



< 



e 



Jr 



[L-^{K)i-Voh, 1)] o ^-J . VUc!vlv,)dSds 

+ CN2{uo,F) + C f K{s)\\W^oh,\\l,,^.ds + {5 + Cf/^) f \Mj,s,n)ds 
Jo Jo 



for some constant C depending on M, S and di. 



52 



C.H. ARTHUR CHENG, DANIEL COUTAND, AND STEVE SHKOLLER 



Appendix C. L^H^ estimates for vt 
By the chain rule and integrating by parts, 



Q[Lj,{K){-Woh,l)]ofi' 



t Jr 



Lhi^T-n) o?7^(-Vo/i0 77^,l) •w^jdS 



+ / Sfjl ■ [Vo[LfXhn)K-Voh, 1)J o fj- ■ v^tdS (= K^) 

e\[LjXhM'^oh,-l)]\ ofj^-v^,dS. i=K2) 
L J t 

The first term is bounded by 

C(M) II S||h3(o) [\\yth.\\mr) + l] l|wKtl|L2(r) • 

After integrating by parts, the most difficult term to estimate in Ki consists of the 
integral 

[Vdet(go)A"'3''^/i«,„^],^5(Vo/i, -1)1 o fj^Vav^^dS. 



/r y/dct{g^ 
Integrating from to t and integrating by parts in time, we find that 

V 



Jr y/det{go) 

t 



[./A^t{^)A''f"''K^^pi^5{Vji,^l) 



/o Jt Vdet(5o) 
where -R3 is bounded by 



[yditM^"^'''/iK,a/3]*,75(Vo/i, -1) 



° ff^oVi^idSds 

o fj'^VoVudSds + i?3 



1 + W^tWHun) \\^oh.\\i2fr)ds + S^WythnWt 



(r) 



+S h>^fHHn)ds + (S + Cf^^) \\v^tfHHn)ds 



C 



for some constant C depending on M, (5 and 52- Next, using that 

[{-Vji, 1) o ,7-] . Vo«« = h\V^Kt) ° V" + h\Vlh o ^^ 0) • v«, 

and integrating by parts, we find that the integral on the right-hand side is identical 
to 



1 



JT 



1 



\jAet{go 



Vdet(go)ev6*A"^^*l h^t^^ph^t,jsdSds + R4 



where 



4j. I|2 J„ I X / ll„. I|2 



|i?4| < CiM)CiS) / ||V^/j«||i.(r)ds + <5 / ||««||^3(o)ds. 
Jo Jo 

By interpolation, the integral part is bounded by 



C 



N{uo,F)+ \\Vthn\\i2^r)ds +S \K\\H3f^,,^ds + Ct ||«.J|,i(o)(^s 



NAVIER-STOKES INTERACTING WITH A FLUID SHELL 

for some constant C depending on M and 5. Therefore, Ki satisfies 



53 



Kids 



<C 



K{s){\\VlK\\h^^) + \\ViK,\\i,^^) 



1 



(is + (52||Vo/i„||^2(r 



(r) 



+ (S + Ct'/^) \KfHS^,,^ds + id + Ct'/^) \\v,t\\hin)ds (C.l) 

Jo Jo 

for some constant C depending on M , S and 52- 

For K2, by time differentiating the evohition equation, we find that 

(-V0/1 o f}'^, l)v^t = h^tt of]^ + v'^ ■ (Vo/iKt) o ^■^ - -D^ • {Vlh o f]^, 0) • v^ 

and hence (after a change of coordinates) 

i^2 = y" [L-^{K)]tKudS + J [LjXKMiv^ o f}-^) ■ {VoKt)]dS (= K3) 

[L-hiKMiVohu 0) • K ° n]dS (= ir4) 



[L-^{KmVoht,0)-{v^t°fj-^)]dS {=Ke). 



r 



For the first term, we have 



J [Lj,{h,)]tKttdS ^\jj A^P-<'K,^^f,K,^^sdS 



h v/det(.go) 
where -R5 is bounded by 



yd^t(^(A"'3^^)J K^^pKudS {=Kr)+R5 (C.2) 



,75 



C 



l + ll^t|lH2.5(r) l + ||V§/i«Ji.(r) +(5 ||««||h2(o) + ||V2z;,||2^i(jj,) 



+ Si\\vKt\\Hi(n) 
for some constant C depending on M, (5 and 61. Also, by the inequahty |j/iKtt|JL4(p) < 
C{M) \\v^\\H^(^Q) + \\vKt\\HUn) , 

1 



\K^\ < C||[v/d^tM(^"^^')*],75lk-- 



= (r) 

iu l|2 



= "-K,Q/3"Ktt 



^0.5(r) 
2 



Remark 21. The bound for K7 can be refined even further as 

\Kr\ < C{M)C{S)\\ht\\l^.r,(r)\\^lh^\\m.^(^^ + S\\v^\\j,s(n) + S\Kt\\m(n); 
it is this inequality that will be used in the proof of the fixed-point argument. 
It remains to estimate K3 to Kq. By proper use of Holder's inequality, 



\K. 



\K. 



iKel < C 






l + llVXlli^rr 



(r) 



+ {s + ct'/^)\K\\U,,. + s\\v^,rH^ 



i(0) 
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for some constant C depending on M and 5. For K^, most of the terms can be 
estimated in the same fashion except the term 



1 



/r \/det(go) 
which is identical to 
1 



^dei{go)A'^^'^'Kt,o.p] [{"^ohtns, 0) • [v^ o rl]dS 



It ^ A/det(go) 
where 



v/d^t(^A"^^*/i«,,„^] [(Vo/i,^5,0) • K ° ^"")]}/5 (= Ks) + 



Rr 



\R&\ < C||/i||^5.5(r) [lkK||i2(j^) + ||Vg/iKJi2(r)J + ^||w«:||h3(o) + SiWv^tWmiQ) 

for some constant C depending on M, S and 6i. Time integrating K^ and use the 
interpolation inequality together with Yomig's inequality, we find that 



Ks{s)ds <C{M) ||uo|Ih2.5(o) + l|Vo^Kt||L2(o)||uK||Li(n) 
< C{M)C{5i,52)N^{u„,F) + 52\\VlKt\\l,^^) + 5i f Wv^tWmmds 



(C.3) 



where 



N3{uo,F) :- ||uo||_H-2E.(o) + l|wo||//i-5(r) + \\P\\L^{o,T-HHn)) 

+ ll-^*llL2(o,T;ffi(f^)') + ll^(0)llffi(f^) + -*- 



and we use ll'yKllfl-i(si) < C 
hence 



\vKt\\Hi{n)ds + lluoWn^n) 



to obtain (|C.3[) . and 



Y. m < c[i + mu.^r) + fhtWUrHT)] [i + ikKiii^(o) + iivi^/^.m^jr)" 

i=3 

+ iS + Ct'/')\Kfm^^^ + SilKtfmm + Ks (C.4) 

with Ks satisfying inequality (|C.3|) . Finally, combining all the estimates, 

[L,- (/i«)(Vo/i, -1)] o f]^] ■ v^tdS + CTVsK, F) 



t rt 

Jo Jr 



+ C Kis) \\v.\\Un) + \\yth.\\i2^r) + mh.t\\i 



Hn 



ds 



(C.5) 



+ (5 + Ct'/') / \\v.\\l.(n)ds + {S,+Ct'/') / \KtrH.(n)d^ + S2\\VtK\ 



for some constant C depending on M, S, Si and 62- 



i(0)^ 



\L^r) 
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